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ABSTRACT 

In  1996,  the  circuit  of  the  low-speed  wind  tunnel  was  modified  by  extending  the  length 
of  the  test  section  and  installing  a  redesigned  contraction.  This  was  done  to  enable 
longer  models  to  be  tested  in  the  tunnel  and  also  to  possibly  obtain  some  improvement 
in  the  quality  of  the  flow  in  the  test  section.  In  this  report,  the  results  are  given  of  a 
detailed  calibration  of  the  flow  in  the  extended  test  section  of  the  tunnel.  The 
calibration  was  performed  to:  (1)  provide  a  data  base  of  the  flow  quality,  to  be  used 
when  analysing  tests  carried  out  in  the  tunnel  since  1996,  and  (2)  provide  a  data  base  for 
assessing  the  changes  in  the  flow  quality  resulting  from  installing  the  extended  test 
section  and  new  contraction  in  1996.  Longitudinal  mean-flow  velocities,  flow  angles, 
turbulence  intensities  and  spectra,  were  measured  at  up  to  155  grid  points  across  the 
flow  for  longitudinal  positions  of  xT  =  -2.0,  -1.0,  0.0  and  +1.0  m  from  the  centre  of  the 
turntables  in  the  test  section  for  nominal  free-stream  velocities  of  30,  45,  60  and  75  m/s. 
The  most  important  area  of  the  test  section,  where  models  are  usually  located,  is  the 
region  comprising  approximately  the  central  50%  of  the  cross  sectional  area  of  the  test 
section.  In  this  region,  over  the  longitudinal  measurement  range  (xT  =  -2.0  to  +1.0  m), 
the  following  flow  non-uniformities  were  observed.  For  the  four  velocities,  longitudinal 
mean-flow  velocities  deviated  within  the  range  -1.2%  to  +1.6%  from  their  values  at  the 
centre  of  the  tunnel  test  section  (xT  =  0.0  m).  For  velocities  of  30  and  60  m/s,  horizontal 
flow  angles  deviated  within  the  range  -0.8°  to  +0.8°  and  vertical  flow  angles  deviated 
within  the  range  -1.2°  to  +0.4°.  For  velocities  of  30  and  60  m/s,  the  u-component 
turbulence  intensities  were  generally  below  0.4%,  and  the  v-  and  w-component 
intensities  were  generally  below  0.7%.  Corresponding  spectra  for  the  u-  v-  and 
w-components  of  the  turbulence  were  relatively  smooth  over  the  frequency  range 
0.3  Hz  to  10  kHz,  with  no  significant  spikes  associated  with  irregularities  in  the  flow. 
There  was  a  slight  improvement  in  the  quality  of  the  flow  in  the  test  section  resulting 
from  the  modifications  to  the  tunnel  in  1996. 
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Calibration  of  the  Flow  in  the  Extended  Test  Section  of  the 
Low-Speed  Wind  Tunnel  at  DSTO 

Executive  Summary 

The  low-speed  wind  tunnel  at  the  Platforms  Sciences  Laboratory  (PSL)  in  the  Defence 
Science  and  Technology  Organisation  (DSTO)  was  commissioned  in  1941.  It  has  a 
closed  return  circuit  with  an  octagonal  test  section  2.74  m  wide  and  2.13  m  high,  and 
the  maximum  velocity  in  the  test  section  is  about  100  m/s.  The  tunnel  circuit  remained 
essentially  unchanged  from  1941  until  1996,  when  the  test  section  length  was  increased 
from  4.12  m  to  6.56  m  and  a  new  contraction  was  installed  upstream  of  the  original 
contraction.  This  was  done  to  enable  longer  models  to  be  tested  and  also  to  possibly 
obtain  some  improvement  in  the  quality  of  the  flow. 

The  flow  in  a  tunnel  test  section  is  never  completely  uniform.  Detailed  measurements 
need  to  be  taken  in  the  empty  test  section  to  see  how  the  flow  deviates  from  location  to 
location.  Tunnels  are  calibrated  to:  (1)  establish  if  corrections  have  to  be  made  to  test 
results,  (2)  indicate  if  adjustments  are  needed  to  the  orientation  of  a  model  when  testing 
and  (3)  assess  the  suitability  of  the  tunnel  for  a  given  test.  Although  some  early 
calibration  data  was  taken  for  specific  purposes,  it  was  not  comprehensive. 

In  this  report,  the  results  are  given  of  a  detailed  and  comprehensive  calibration  of  the 
flow  in  the  extended  test  section  of  the  tunnel.  The  calibration  was  performed  to: 

(1)  provide  a  data  base  to  facilitate  detailed  analysis  of  tests  conducted  since  1996,  and 

(2)  provide  a  data  base  for  assessing  the  flow  quality  resulting  from  installing  an 
extended  test  section  and  a  new  contraction  in  1996.  A  range  of  longitudinal  mean-flow 
velocities,  flow  angles,  turbulence  intensities  and  spectra  were  measured  in  the  test 
section.  Data  were  taken  at  up  to  155  grid  points  across  the  flow  at  each  of  four 
longitudinal  distances  of  xT  =  -2.0,  -1.0,  0.0  and  +1.0  m  from  the  centres  of  the  upper 
and  lower  turntables  and  for  nominal  free-stream  velocities  of  30, 45,  60  and  75  m/s. 

The  most  important  area  of  the  test  section,  where  models  are  usually  located,  is 
approximately  the  central  0.7  of  the  width  and  0.7  of  the  height,  i.e.  the  central  50%  of 
the  cross  sectional  area  of  the  test  section.  In  this  region,  over  the  longitudinal  measure¬ 
ment  range  xT  =  -2.0  to  +1.0  m,  the  following  flow  non-uniformities  were  observed. 
For  velocities  of  30,  45,  60  and  75  m/s,  longitudinal  mean-flow  velocities  deviated 
within  the  range  -1.2%  to  +1.6%  from  their  values  at  the  centre  of  the  tunnel  test 
section  (xT  =  0.0  m).  For  velocities  of  30  and  60  m/s,  horizontal  flow  angles  deviated 
within  the  range  -0.8°  to  +0.8°  and  vertical  flow  angles  deviated  within  the  range  -1.2° 
to  +0.4°.  These  are  the  maximum  limits  of  variation  over  the  complete  velocity  range, 
but  for  any  given  velocity  the  deviations  in  mean  velocities  and  flow  angles  were 
generally  substantially  less  than  these  values.  For  velocities  of  30  and  60  m/s,  the 
u-component  turbulence  intensities  were  below  0.4%,  and  the  v-  and  w-component 
intensities  were  below  0.7%.  Spectra  for  the  u-  v-  and  w-components  of  the  turbulence 
were  relatively  smooth  over  the  frequency  range  0.3  Hz  to  10  kHz,  with  no  significant 
spikes  associated  with  irregularities  in  the  flow.  The  quality  of  the  flow  in  the  extended 
test  section  was  found  to  be  slightly  better  than  that  in  the  original  test  section. 

The  literature  indicates  that,  for  a  low-speed  tunnel,  longitudinal  velocities  should  not 
deviate  from  the  mean  value  by  more  than  ±0.75%,  and  horizontal  and  vertical  flow 
angles  should  not  deviate  from  the  longitudinal  direction  by  more  than  ±0.25°.  The 
u-component  turbulence  intensities  should  be  less  than  about  0.1%  for  research,  but 
intensities  as  high  as  0.5%  are  acceptable  for  developmental  testing.  The  variations  in 
mean  velocities  and  flow  angles  for  the  DSTO  tunnel  are  a  little  higher  than  some  other 
tunnels.  In  the  future,  additional  modifications  could  be  made,  such  as  replacing  the 
honeycomb  or  installing  screens,  to  improve  the  quality  of  the  flow  in  the  test  section. 
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Notation 

Cp  Pressure  coefficients  associated  with  the  six-hole  probe, 

Cpn  =  (  Pn  -  Pref  )7(0-5  P  Vy),  (n  =  1  to  6). 

f  Frequency,  (Hz). 

H  Height  of  the  test  section  of  the  low-speed  wind  tunnel,  (H  =  21 34  mm). 

M  Molecular  weight  of  air  (M  =  28.9644  ),  (dimensionless). 

pcr,  p^  Static  pressures  at  the  upstream  and  downstream  ends  respectively  of  the 
tunnel  contraction,  (Pa). 

pref  Reference  pressure  (Pa). 

pp  p„  p, ...  Pressures  associated  with  the  six-hole  probe  or  the  cross-tunnel  rake 

(pressures  are  not  necessarily  total  or  static),  (Pa). 

R*  Universal  gas  constant  (R*  =  8314.32  J/[K.kg  Mol]). 

Su,  Sv,  Sw  Normalized  power  spectral  densities  for  the  u,  v  and  w  components  of  the 
turbulence  respectively,  (dimensionless). 

t  Static  temperature  of  the  air  in  the  test  section,  (K). 

u  Fluctuating  component  of  velocity  in  the  xT  direction,  (m/s). 

U  Velocity  in  the  xT  direction  (U  =  U  +  u),  (m/s). 

U  Mean  velocity  in  the  xT  direction,  (m/s). 

UB  Mean  velocity  in  the  xB  direction  in  the  probe  body  coordinate  system,  (m/s) 

UNOM  Nominal  free-stream  velocity  in  the  xT  direction,  (m/s). 

UREF  Reference  free-stream  velocity  at  the  “centre”  of  the  test  section,  (m/s). 

UT  Mean  velocity  in  the  xT  direction  in  the  tunnel  coordinate  system,  (m/s). 

U_,0  Mean  velocity  in  the  xT  direction  at  the  location  where  xT  =  -2.0  m, 

yT  =  0.0  m  and  zT  =  0.0  m,  (m/s). 

U_10  Mean  velocity  in  the  xT  direction  at  the  location  where  xT  =  -1.0  m, 
yT  =  0.0  m  and  zT  =  0.0  m,  (m/s). 

Unn  Mean  velocity  in  the  xT  direction  at  the  location  where  xT  =  0.0  m, 

yT  =  0.0  m  and  zT  =  0.0  m,  (m/s). 

U,n  Mean  velocity  in  the  xT  direction  at  the  location  where  xT  =  1.0  m, 

yT  =  0.0  m  and  zT  =  0.0  m,  (m/s). 

v  Fluctuating  component  of  velocity  in  the  yT  direction,  (m/s). 

VB  Mean  velocity  in  the  yB  direction  in  the  probe  body  coordinate  system,  (m/s) 

VT  Mean  velocity  in  the  yT  direction  in  the  tunnel  coordinate  system,  (m/s). 
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VY  Velocity  vector  as  measured  with  the  six-hole  probe  (yawmeter) 

(Figure  14),  (m/s). 

w  Fluctuating  component  of  velocity  in  the  2^  direction,  (m/s). 

W  Width  of  the  test  section  of  the  low-speed  wind  tunnel,  (W  =  2743  mm). 

WB  Mean  velocity  in  the  z„  direction  in  the  probe  body  coordinate  system,  (m/s) 

WT  Mean  velocity  in  the  zT  direction  in  the  tunnel  coordinate  system,  (m/s). 

xu  yB  z«  B°dy  coordinate  system  (right-handed)  for  the  six-hole  probe  (Figure  14). 

The  origin  is  located  at  centre  of  the  hemisphere  which  forms  the  tip  of  the 
probe.  The  xH  axis  is  along  the  probe  longitudinal  axis,  the  yn  axis  is  in  the 
plane  passing  through  holes  2,  4  and  5  and  the  z„  axis  is  in  the  plane 
passing  through  holes  1,  3  and  5.  The  axes  remain  fixed  with  respect  to 
the  probe. 

V  yB»  zi.  Directions  or  distances  in  the  probe  body  coordinate  system  (Figure  14),  (m) 

xT  yT  zt  Tunnel  coordinate  system  (right-handed)  (Figure  2). 

The  origin  is  located  at  the  midpoint  of  the  line  joining  the  centres  of  the 
lower  and  upper  turntables  in  the  test  section.  The  xT  axis  is  on  the  test 
section  longitudinal  centreline  and  is  positive  in  the  downstream  direction, 
the  yT  axis  is  horizontal  and  is  positive  to  port  and  the  z^  axis  is  positive 
vertically  downwards.  The  axes  remain  fixed  with  respect  to  the  tunnel. 

xT»  y-p  zT  Directions  or  distances  in  the  tunnel  coordinate  system  (Figure  2),  (m). 


Horizontal  flow  angle  in  the  tunnel  coordinate  system  (Figure  19),  (deg). 

Vertical  flow  angle  in  the  tunnel  coordinate  system  (Figure  19),  (deg). 

Pitch  angle  for  the  six-hole  probe,  relative  to  the  tunnel  coordinate  system 
x^^  (Figure  14),  (deg). 

Flow  angle  as  measured  with  the  six-hole  probe  (yawmeter)  (Figure  14) 
(deg).  '  '  “ 

Matrix  of  direction  cosines  used  to  transform  velocities  from  the  tunnel 
coordinate  system  to  the  probe  body  coordinate  system. 

Matrix  of  direction  cosines  used  to  transform  velocities  from  the  probe  body 
coordinate  system  to  the  tunnel  coordinate  system. 

Density  of  the  air  in  the  test  section,  (kg/m3). 

Roll  angle  for  the  six-hole  probe,  relative  to  the  tunnel  coordinate  system, 
xTyTZr  (Figure  14),  (deg). 

Flow  angle  as  measured  with  the  six-hole  probe  (yawmeter)  (Figure  14) 
(deg).  '  " 

Yaw  angle  for  the  six-hole  probe,  relative  to  the  tunnel  coordinate  system 
(Figure  14),  (deg). 
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1.  Introduction 

The  low-speed  wind  tunnel  (LSWT)  at  the  Platforms  Sciences  Laboratory  (PSL)  in  the 
Defence  Science  and  Technology  Organisation  (DSTO)  was  designed  in  1939  and 
commissioned  in  1941  (Reference  1).  It  has  a  closed  return  circuit  with  an  octagonal 
test  section  2.74  m  wide  and  2.13  m  high,  and  the  maximum  velocity  obtainable  in  the 
test  section  is  about  100  m/s.  At  the  downstream  end  of  the  test  section  there  is  a 
pressure  equalizing  slot  which  maintains  the  test  section  at  near  atmospheric  pressure. 
The  flow  circuit  remained  essentially  unchanged  from  1941  until  1996,  when  the  length 
of  the  test  section  was  increased,  enabling  longer  models  to  be  tested,  and  a  redesigned 
contraction  was  installed.  Within  the  constraints  imposed  by  the  space  available  in  the 
existing  settling  chamber  and  the  new  contraction  design,  the  length  of  the  test  section 
was  increased  from  4.12  m  to  6.56  m.  The  modifications  to  the  tunnel  in  1996  followed 
the  design  recommendations  given  by  Watmuff  (Reference  2).  The  modified  tunnel  is 
shown  diagrammatically  in  Figure  1  (Reference  1),  which  also  shows  an  outline  of  the 
original  contraction.  Details  of  the  new  contraction  design  are  given  in  Appendix  A. 

Ideally,  the  flow  in  the  test  section  of  a  wind  tunnel  is  uniform  throughout,  is  parallel  to 
the  tunnel  walls  (for  parallel-walled  tunnels)  and  has  zero  turbulence.  However,  the 
flow  is  never  perfect  and  it  is  necessary  to  calibrate  the  flow  to  determine  the  extent  of 
the  irregularities.  To  calibrate  the  flow  in  a  tunnel,  detailed  flow  measurements  are 
taken  throughout  the  empty  test  section  to  see  how  the  flow  deviates  from  location  to 
location.  The  calibration  establishes  the  corrections  to  be  made  to  test  results,  e.g. 
corrections  may  have  to  be  made  to  measurements  if  the  angle  of  the  flow  incident  on 
one  wing  of  an  aircraft  is  slightly  different  from  that  on  the  other  wing.  The  calibration 
also  indicates  what  adjustments  have  to  be  made  to  the  orientation  of  a  model  when 
testing,  e.g.  the  flow  may  have  some  angularity  so  that  a  strain-gauge  balance  in  a 
model  may  need  to  be  aligned  with  the  flow  rather  than  the  axes  of  the  tunnel.  A 
calibration  enables  an  assessment  to  be  made  of  the  suitability  of  a  tunnel  for  a  given 
project,  e.g.  a  tunnel  having  a  high  free-stream  turbulence  intensity  may  be  unsuitable 
for  fundamental  experiments  in  which  boundary-layer  transition  is  important,  since  high 
turbulence  levels  can  significantly  affect  the  characteristics  of  the  flow  being  studied.  A 
calibration  establishes  the  limitations  of  the  tunnel. 

A  comprehensive  calibration  of  the  flow  in  the  test  section  of  the  tunnel  with  the 
extended  test  section  and  new  contraction  has  been  carried  out  and  details  of  the 
calibration  are  given  in  this  report.  It  is  important  that  the  tunnel  is  calibrated  to  ensure 
that  the  quality  of  the  data  produced  is  suitable  for  its  intended  purpose.  An  extensive 
range  of  longitudinal  mean-flow  velocities,  flow  angles,  turbulence  intensities  and 
turbulence  spectra  were  measured.  Data  were  taken  at  up  to  155  grid  points  across  the 
flow  at  each  of  four  longitudinal  distances  of -2.0,  -1.0,  0.0  and  +1.0  m  from  the  datum 
of  the  tunnel  coordinate  system  (see  Section  2.1)  and  for  nominal  free-stream  velocities 
of  30,  45,  60  and  75  m/s.  The  calibration  data  are  presented  and  analysed  and 
explanations  of  flow  behaviour  are  given. 
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An  earlier  calibration  was  carried  out  not  long  before  the  extended  test  section  and  new 
contraction  were  installed  in  1996.  The  calibration  measurements  taken  were  similar  to 
those  listed  above.  Details  of  this  earlier  calibration  are  given  in  Reference  3.  Selected 
data  from  the  earlier  calibration  are  given  in  the  current  report  to  assess  the 
effectiveness  of  the  extended  test  section  and  the  new  contraction.  All  comparisons  will 
be  made  for  the  flow  cross  section  at  the  datum  of  the  tunnel  coordinate  system  (see 
Section  2.1). 


2.  Mean- Velocity  Measurements 


2.1  Test  Schedule 

Longitudinal  mean  velocities,  U,  were  measured  in  the  test  section  of  the  tunnel  at  the 
155  grid  points  shown  in  Figure  2  for  each  of  four  longitudinal  positions, jcT  (see 
below),  of  -2.0,  -1.0,  0.0  and  +1.0  m  and  for  nominal  free-stream  velocities,  UN0M,  of 
30, 45,  60  and  75  m/s. 

The  origin  of  the  tunnel  coordinate  system,  xT  yTzT,  shown  in  Figure  2,  is  located  at  the 
midpoint  of  the  line  joining  the  centres  of  the  concentric  lower  and  upper  turntables  in 
the  test  section.  The  xT  axis  is  on  the  test  section  longitudinal  centreline  and  is  positive 
in  the  downstream  direction,  the  yT  axis  is  horizontal  and  is  positive  to  the  port  side  of 
the  tunnel  and  the  zT  axis  is  positive  vertically  downwards.  The  axes  form  a  right- 
handed  coordinate  system  and  remain  fixed  with  respect  to  the  tunnel. 

In  Figure  2,  the  grid  points  are  spaced  at  intervals  of  137.0  mm  in  the  yT  direction, 
corresponding  to  the  spacing  of  the  pressure  probes  on  the  cross-tunnel  rake 
(Section  2.3),  and  at  intervals  of  213.4  mm  in  the  zT  direction.  The  yT  spacing  corres¬ 
ponds  to  l/20th  of  the  width  of  the  octagonal  test  section  and  the  zT  spacing  corresponds 
to  l/10th  of  the  height  of  the  test  section  (W  =  2743  mm  and  H  =  2134  mm). 


2.2  Procedure  Used  to  Measure  Reference  Free-Stream  Velocities 

Reference  free-stream  velocities  at  the  “centre”  of  the  test  section  of  the  tunnel  were  set 
using  the  static  pressures,  pcu  and  pCD,  at  the  upstream  and  downstream  ends  respectively 
of  the  new  contraction.  The  new  contraction  was  calibrated  shortly  after  being  installed 
and  it  was  determined  that 


Urrf 


1,079  x2(pcu-pCD) 
P 


(1) 
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Figure  2.  Grid  points,  shown  circled,  at  which  mean  velocities  and  flow  angles  were  measured. 


The  1.079  factor  is  an  empirical  calibration  factor  that  applies  specifically  to  the  LSWT. 
The  density  of  the  air  in  the  test  section,  p,  is  given  by  (Reference  4) 


(2) 


M  is  the  molecular  weight  of  air  (M  =  28.9644  [dimensionless]),  R*  is  the  universal  gas 
constant  (R*  =  8314.32  J/[K.kg  Mol])  and  t  is  the  static  temperature  of  the  air  in  the  test 
section.  Equation  1  was  used  to  set  the  reference  free-stream  velocity  in  the  test  section 
when  the  detailed  calibration  measurements  were  taken. 
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2.3  Cross-Tunnel  Rake 

The  cross-tunnel  rake  shown  in  Figure  3  was  used  to  measure  pressures  at  the  different 
grid  points  shown  in  Figure  2.  The  rake  is  made  up  of  a  number  of  sections  bolted 
together  and  contains  15  total-pressure  probes  and  4  static-pressure  probes,  arranged  as 
shown  in  Figure  3.  For  z/H  locations  of  0.0,  ±0.1  and  ±0.2  (Figure  2),  the  rake 
contained  probes  1  to  19  and  was  bolted  to  the  side  walls  of  the  test  section  using 
mounting  brackets.  To  allow  for  the  corner  fillets  in  the  octagonal  test  section,  the 
length  of  the  rake  could  be  altered  by  using  different  end  sections.  For  zr/H  =  ±0.3,  the 
rake  contained  probes  2  to  18  and  for  z/H  =  ±0.4,  the  rake  contained  probes  4  to  16. 
When  the  rake  was  shortened,  supporting  stands  were  attached  to  the  rake  so  that  it 
could  be  bolted  to  the  floor  or  the  roof  of  the  test  section. 

2.4  Procedure  Used  to  Measure  Mean  Velocities 

Longitudinal  mean  velocities  at  different  grid  points  were  measured  using  the  pressure 
probes  on  the  cross-tunnel  rake.  Pressures  were  measured  using  a  Pressure  Systems 
Incorporated  (PSI)1  8400  Electronic  Measurement  System,  configured  with  a  32  port 
Electronic  Pressure  Scanner  capable  of  measuring  32  differential  pressures.  The 
scanner  had  a  range  of  ±2490  Pa  (±10.0  inches  of  water)  for  UN0M  =  30  m/s  and  a  range 
of  ±6895  Pa  (±1.0  psi)  for  UNOM  =  45,  60  and  75  m/s.  Software  used  to  control  and 
operate  the  system  is  located  on  a  PC.  The  features  of  the  system  are  described  in 
References  5  and  6. 

The  total  and  static  pressure  ports  on  the  probes  on  the  rake,  as  well  as  the  static 
pressure  ports  in  the  upstream  and  downstream  ends  of  the  contraction,  used  when 
measuring  pcu  and  pCD  (static  pressure  in  the  test  section)  respectively,  were  connected 
to  ports  on  the  scanner,  and  all  ports  were  referenced  to  pCD.  Prior  to  any  given  run,  the 
pressure  transducers  in  the  scanner  were  calibrated  over  the  range  of  pressures  likely  to 
be  encountered  during  subsequent  measurements.  Four  different  ranges  of  pressures 
were  used,  corresponding  to  nominal  free-stream  velocities  of  30,  45,  60  and  75  m/s. 
The  accuracy  of  the  calibration  was  verified  by  checking  that  all  measured  differential 
pressures  were  close  to  zero,  as  observed  graphically  on  a  monitor,  when  there  was  no 
flow  in  the  tunnel. 

The  software  allowed  the  user  to  specify  sampling  controls  when  acquiring  data.  The 
differential  pressures  corresponding  to  each  of  the  ports  were  sampled  50  times  at 
intervals  of  65  000  (is  and  then  average  pressures  for  each  port  were  calculated  and 
recorded.  The  65  000  ps  interval  was  the  maximum  allowable  with  the  system  and 
corresponds  to  a  sampling  frequency  of  15.4  Hz.  The  total  sampling  time  for  the  50 
samples  was  about  3.2  s.  By  setting  the  sampling  frequency  at  the  lowest  possible  value 
allowed  by  the  equipment,  any  large  scale,  slowly  varying,  pressure  oscillations  in  the 
tunnel  were  more  likely  to  be  averaged.  The  above  sampling  was  repeated  10  times 

1  Pressure  Systems  Incorporated,  34  Research  Drive,  Hampton,  VA,  23666,  USA. 
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4.93 


Figure  3b.  Cross-tunnel  rake  used  to  measure  mean-flow  velocities  showing  diagrammatic 
representation  of  rake  (top),  cross  section  of  rake  (middle)  and  details  of  probes  (bottom). 
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with  a  5  s  time  delay  between  sampling  each  data  set.  The  differential  pressures  corres¬ 
ponding  to  each  of  the  ports  were  then  averaged  for  the  10  sets  of  data  to  obtain  overall 
averages,  each  of  which  corresponds  to  50  x  10  =  500  individual  samples. 

As  can  be  seen  from  Figure  3,  the  arrangement  of  the  probes  on  the  rake  meant  that 
either  the  total  pressure  or  the  static  pressure  was  measured  at  each  grid  point.  Since  it 
was  necessary  to  know  both  the  static  and  the  total  pressure  at  a  given  point  to  calculate 
the  velocity  at  that  point,  approximations  had  to  be  made  when  calculating  velocities. 
When  calculating  velocities  at  grid  points  where  total  pressures  were  known,  nearby 
static  pressures  were  used.  For  example,  the  velocity  at  point  4  was  based  upon  p4  (total 
pressure),  p,  (static  pressure)  and  p7  (static  pressure).  A  mean  velocity  denoted  by  U4A 
was  computed  using 


U 


4A 


2(p4  ~  pJT 3 
p 


(3) 


and  a  mean  velocity  denoted  by  U4lt  was  computed  using 


U 


4H 


2(p4-p7)l0-5 

P 


(4) 


The  mean  velocity  U4  at  position  4  was  then  determined  from  U4A  and  U4B  using  linear 
interpolation  as  follows 


U.  =  0.75  [UjJ  +  0.25  [U4II] 


(5) 


At  grid  points  where  static  pressures  were  known,  mean  velocities  were  determined 
using  linear  interpolation  of  the  calculated  adjoining  velocities.  For  example,  the  mean 
velocity  at  point  3  was  determined  from  U2  and  U4  using  linear  interpolation  as  follows 

U,  =  0.5  [U2  +  Uj  (6) 


The  reference  free-stream  velocity,  URI3„  determined  using  equation  1,  varied  slightly 
for  the  different  positions  of  the  cross-tunnel  rake,  due  to  the  practical  difficulties  of  a 
tunnel  operator  setting  the  free-stream  velocity  at  precisely  the_same  value  for  the 
different  rake  positions.  To  obtain  a  consistent  set  of  values  of  U  measured  with  the 
rake,  the  values  of  U  were  non-dimensionalised  with  respect  to  the  appropriate  value  of 
UKI;, .  To  obtain  a  consistent  set  of  results  with  respect  to  a  given  cross  section  in  the 
tunnel,  for  example  at  the  centre  oHhe_cross  section  where  xT  =  1 .0,  the  results  were 
then  normalized  with  respect  to  the  U,  n  /URI;):  ratio  at  this  point. 


8 


DSTO-TR-1384 


2.5  Analysis  of  Mean  Velocities 

Longitudinal  mean  velocities  are_presented  as  contours  of  (U  -  U_20)/U_20, 
(U-U_J/U_10,  (U  -  U00)/U0  „  _and  (U -Uj/U10_(expressed  as  a  percentage)  for 
different  values  of  xT  and  UN0M.  U_20,  U_,  0,  U0 0  and  U,  0  are  the  mean  velocities  in  the  xT 
direction  on  the  tunnel  longitudinal  centreline  (yT  =  0.0  m  and  zT  =  0.0  m)  at  the 
locations  where  xT  =  -2.0,  -1.0,  0.0  and  +1.0  m  respectively.  Contour  plots  of  mean 
velocities,  as  well  as  contour  plots  of  flow  angles  and  turbulence  intensities,  are 
presented  with  the  direction  of  the  flow  out  of  the  page. 

The  most  important  area  of  the  test  section,  where  models  are  usually  located,  is  the 
region  comprising  approximately  the  central  0.7  of  the  width  and  0.7  of  the  height,  i.e. 
approximately  the  central  50%  of  the  cross-sectional  area  of  the  test  section.  Discussion 
on  the  variation  of  mean  velocities,  as  well  as  on  flow  angles  and  turbulence  intensities, 
will  mainly  apply  to  these  central  regions  of  flow  cross  sections. 

2.5. 1  Analysis  With  Respect  to  Each  Cross  Section 

Contours  of  (U  -  U_20)/U_20  for  _xT  =_  -2d)  m,  (U  -  U_10)/U_L0  for  x^.  =  —1.0  m, 
(U-Uj/U0,  for  xT  =  0.0  m  and  (U  -  U,  J/U, „  for  xT  =  +1.0  m,  all  for  UN0M  =  30  m/s, 
are  shown  in  Figure  4.  Corresponding  plots  for  UN0M  =  45,  60  and  75  m/s  are  shown  in 
Figures  5,  6  and  7  respectively.  For  each  of  the  16  sets  of  contours  shown  in  Figures  4 
to  7,  the  normalizing  velocity,  i.e.  U_20>  U_10,  U00  or  UJ  0,  is  associated  with  the  xT 
location  of  the  contours,  i.e.  U_20  is  used  for  contours  at  xT  =  -2.0  m,  and  likewise  for 
U_10?  Uo  n  and  U10.  For  each  of  the  four  values  of  UNOM,  the  general  pattern  of  the 
contours  at  the  four  xT  locations  is  similar,  with  little  evidence  of  the  flow  becoming 
more  uniform  as  it  moves  downstream.  Regions  of  high  or  low  mean  velocity  are 
simply  convected  downstream  with  little  transverse  movement,  indicating  that  there  are 
no  large-scale  swirls  in  the  test  section.  For  UNOM  =  30  m/s,  for  xT  =  -2.0,  -1.0,  0.0  and 
+  1.0  m,  longitudinal  mean-flow  velocities  typically  deviated  between  -0.8%  and  +1.2% 
from  their  values  at  the  centre  of  each  tunnel  cross  section  (Figure  4).  Corresponding 
deviations  forUNOM  =  45,  60  and  75  m/s  were  from  -0.6%  to  +1.2%  (Figure  5),  from 
-0.4%  to  +1.2%  (Figure  6)  and  from  -0.4%  to  +1.4%  (Figure  7)  respectively. 

2.5.2  Analysis  With  Respect  to  the  Central  Cross  Section 

Contours  of  (U  -  U00)/U00  for  UN0M  =  30  m/s  and  xT  =  -2.0,  -1.0,  0.0  and  +1.0  m  are 
shown  in  Figure  8.  Corresponding  plots  for  UN0M  =  45,  60  and  75  m/s  are  shown  in 
Figures  9,  10  and  1 1  respectively.  For  the  16  sets  of  contours  shown  in  Figures  8  to  1 1, 
a  single  normalizing  velocity  at  the  centre  of  the  test  section,  i.e.  U00,  is  used  for  each 
value  of  UN0M,  which  means  that  the  normalizing  velocity  is  not  always  associated  with 
the  xT  location  of  the  contours,  e.g.  U00  is  used  for  contours  at  xT  =  -2.0,  -1.0  and 
+1.0  m.  The  contours  in  each  of  Figures  8  to  11  therefore  show  how  mean  velocities 
vary  in  the  longitudinal  or  xT  direction.  Mean  velocities  at  given  values  of  yT  and  zT 
increase  as  the  flow  moves  downstream,  which  is  consistent  with  the  fact  that  the  test 
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Figure  5  cont'd.  Contours  of  velocity  deviation  for  UN0M  =  45  m/s  for  xT  =  -2.0,  -1 .0, 0.0  and  1.0  m. 

(b)  xT=  0.0  and  1.0  m. 

The  flow  is  out  of  the  page. 
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(b)  xT=  0.0  and  1 .0  m. 


The  flow  is  out  of  the  page. 
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Figure  7  cont'd.  Contours  of  velocity  deviation  for  UN0M  =  75  m/s  for  xT  =  -2.0,  -1 .0, 0.0  and  1 .0  m. 

(b)  xT=  0.0  and  1 .0  m. 

The  flow  is  out  of  the  page. 


Figure  9  cont'd.  Contours  of  velocity  deviation  for  UN0M  =  45  m/s  for  xT  =  -2.0,  -1 .0,  0.0  and  1 .0  m. 

(b)  xT=  0.0  and  1 .0  m. 

The  flow  is  out  of  the  page. 


Figure  10  cont’d.  Contours  of  velocity  deviation  for  UNOM  =  60  m/s  for  xT  =  -2.0,  -1 .0, 0.0  and  1 .0  m. 

(b)  xT=  0.0  and  1 .0  m. 

The  flow  is  out  of  the  page. 


Figure  1 1  cont'd.  Contours  of  velocity  deviation  for  UNOM  =  75  m/s  for  xT  =  -2.0,  -1 .0,  0.0  and  1 .0  m. 

(b)  xT=  0.0  and  1 .0  m. 

The  flow  is  out  of  the  page. 
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section  has  a  fixed  cross-sectional  area  in  the  direction  of  the  flow  so  that  the  effective 
flow  area  reduces  withjncreasing  xT  as  the  wall  boundary  layer  thickens  after  the  end  of 
the  contraction.  ForUNOM  =  30  m/s,  between  xT  =  -2.0  and  +1.0  m,  longitudinal  mean 
velocities  deviated  between  —1.0%  and  +1.4%  from  the  value  at  the  centre  of  the  test 
section  (xT  —  0.0^  yT  =  0.0  and  zT  =  0.0  m)  (Figure  8).  Corresponding  deviations  for 
other  values  of  UN0M  were  from  -1.2%  to  +1.6%  for  UNOM  =  45  m/s  (Figure  9),  from 
I* -°%  t0  +1-4%  for  UN0M  =  60  m/s  (Figure  10)  and  from  -1.2%  to  +1.4%  for 
UNOM  = 75  m/s  (Figure  11).  However,  between  xT  =  -1.0  and  +1.0  m^he  deviations 
were  significantly  smallerjmd  they  ranged  from  -1.0%  to  +1.4%  for  UNOM  =  30  m/s, 
from  -0.6%  to  +1.6%  for_UNOM  =  45  m/s,  from  -0.6%  to  +1.4%  for  U  *=60  m/s  and 
from  -0.8%  to  + 1 .4%  for  UNOM  =  75  m/s. 

2.5.3  Longitudinal  Velocity  Variation 

Selected  data  associated  with  the  contours  shown  in  Figures  8  to  1 1  have  been  replotted 
in  Figure  12  to  indicate  how  the  mean  velocities  in  the  vertical  and  horizontal  longi¬ 
tudinal  centreplanes  vary  in  the  longitudinal  or  xT  direction.  Data  corresponding  to  the 
xT-Zr  plane  that  passes  through  yT  =  0.0  m  (i.e.  the  vertical  plane  that  is  aligned  with  and 
passes  through  the  tunnel  longitudinal  axis)  are  shown  for  z/H  =  0.0,  ±0.2  and  ±0.4. 
Similarly,  data  corresponding  to  the  xT-yT  plane  that  passes  through  zT  =  0.0  m  (i.e.  the 
horizontal  plane  that  is  aligned  with  and  passes  through  the  tunnel  longitudinal  axis)  are 
shown  for  y^W  =  0.0,  ±0.15,  ±0.30  and  ±0.45.  In  Figure  12,  data  corresponding  to  the 
longitudinal  centreline  of  the  tunnel  are  represented  by  thick  lines.  The  longitudinal 
mean  velocities  increase  slightly  in  the  xT  direction,  for  reasons  given  above.  Plots  of 
mean  velocity  for  the  xT-yT  plane  show  that  velocities  on  the  inner  side  of  the  tunnel 
flow  circuit,  i.e.  at  y/W  =  -0.45,  are  the  smallest  and  velocities  on  the  outer  side  of  the 
circuit,  i.e.  at  y/W  =  0.45,  are  the  largest,  for  all  four  values  of  UNOM,  as  would  be 
expected.  As  expected,  mean  velocities  in  the  xT-z.r  plane  do  not  exhibit  any  asymmetry 
for  z/H  =  -0.4  to  0.4,  for  all  four  values  of 

NUM 

2.5.4  Effect  of  Change  in  Nominal  Free-Stream  Velocity 

The  effect  of  change  in_nominaljfree-stream  velocity  are  shown  in  Figures  8  to  1 1.  For 
example,  contours  of  (U  -  U  J/Uao  for  xT  =  0.0  m  and  UNOM  =  30,  45,  60  and  75  m/s  are 
shown  by  the  third  contour  plot  in  each  of  Figures  8  to  11.  Such  a  set  of  contour  plots 
indicates  how  the  flow  at  a  fixed  xT  location  changes  as  UN0M  increases.  The  pattern  of 
the  contours  and  the  magnitudes  of  the  velocity  deviations  are  approximately  the  same 
in  each  case,  indicating  that  the  distribution  of  the  flow  variations  does  not  change 
significantly  with  UN0M. 

2.5.5  Effect  of  Probe  Misalignment  With  Respect  to  the  Flow 

The  probability  that  the  flow  in  the  tunnel  may  be  at  a  slight  angle  to  the  longitudinal 
axes  of  the  probes  on  the  cross-tunnel  rake  (Section  2.3)  will  not  affect  the  test  results 
significantly.  In  Reference  7  it  is  shown  that,  for  total -pressure  and  static-pressure 
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probes  similar  to  those  used  on  the  cross-tunnel  rake  (see  Figure  3),  errors  in  measured 
total  pressures  and  static  pressures  due  to  misalignment  of  the  probes  with  the  flow  are 
negligible  for  flows  whose  angle  of  incidence  is  within  ±5°  of  the  longitudinal  axis  of  a 
total-pressure  probe  and  within  ±1°  of  the  longitudinal  axis  of  a  static-pressure  probe. 
The  angularity  of  the  flow  in  the  test  section  is  generally  less  than  1°  (see  Section  3.7), 
so  that  errors  in  the  dynamic  pressure  and  hence  the  velocity  will  be  small,  and  although 
strictly  speaking  the  contours  shown  in  Figures  4  to  11  depict  slightly  inclined 
velocities,  they  can  be  considered  to  represent  longitudinal  velocities. 

General  discussion  on  the  possible  causes  of  the  flow  non-uniformities  in  the  test 
section  of  the  LSWT  and  suggested  ways  of  improving  the  flow  are  given  in  Section  8. 

2.6  Errors  in  Measured  Mean  Velocities 

There  are  small  errors  in  the  velocity  deviations,  such  as  (U-U00)/Uan,  due  to 
instrumentation  errors.  It  was  necessary  to  measure  air  temperatures  and  differential  air 
pressures  when  determining  mean  velocities.  Errors  in  the  measured  temperatures  did 
not  affect  the  velocity  deviations  as  they  are  expressed  as  non-dimensional  velocity 
ratios.  The  indicated  readings  of  the  scanners  used  to  measure  the  differential  pressures 
were  accurate  to  within  ±0.05%  of  their  full-scale  reading.  Three  different  types  of 
scanners  were  used  when  measuring  pressures  and  their  ranges  were  ±2490  Pa  (±10.0 
inches  of  water),  ±6895  Pa  (±1.0  psi)  and  159  kPa  (23X)  psi).  For  UNOM  =  30,  45,  60  and 
75  m/s,  errors  in  velocity  deviations  (e.g.  (U-U00)/U00)  were  about  ±0.2,  ±0.3,  ±0.2 
and  ±0. 1  percentage  points  respectively. 


3.  Flow- Angle  Measurements 


3.1  Test  Schedule 

Flow  angles  were  measured  in  the  test  section  at  the  155  grid  points  shown  in  Figure  2 
for  xT  =  -2.0,  -1.0,  0.0  and  +1 .0  m  and  for  UNOM  =  30  and  60  m/s. 


3.2  Six-Hole  Probe 

A  six-hole  hemispherical  head  yawmeter,  designed  and  calibrated  by  Melbourne  & 
Griss  (Reference  8),  was  used  to  measure  magnitudes  of  mean  velocities  and 
corresponding  flow  angles  at  the  different  grid  points.  The  probe  is  depicted  diagram- 
matically  in  Figure  13.  Holes  1  to  5  are  located  in  the  hemispherical  head  of  the  probe 
as  shown  and  the  series  of  20  interconnected  holes,  designated  hole  6,  are  equally 
spaced  around  the  circumference  of  the  probe  and  are  located  22.2  mm  (3.5  probe 
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Figure  13.  Six-hole  probe  (based  on  a  diagram  given  in  Reference  8). 


diameters)  from  the  probe  tip.  Each  of  holes  1  to  6  is  connected  internally  to  ports  at 
the  rear  of  the  probe  to  enable  the  individual  pressures  to  be  measured.  A  datum  line 
scribed  on  the  probe  body  near  hole  1  and  contained  in  the  plane  passing  through  holes 
1,  3  and  5,  was  used  to  differentiate  the  holes  and  to  set  the  probe  roll  angle  in  the 
tunnel  (Section  3.4). 

The  direction  of  the  velocity  vector,  Vv,  relative  to  the  probe  body  coordinate  system, 
x»yHZB, is  defined  in  terms  of  flow  angles  0Y  and  c(>v,  as  shown  in  Figure  14.  The  coord¬ 
inate  system  remains  fixed  with  respect  to  the  probe  and  its  origin  is  located  at  the 
centre  of  the  hemisphere  which  forms  the  tip  of  the  probe.  The  x15  axis  is  along  the 
probe  longitudinal  axis,  the  ylt  axis  is  in  the  plane  passing  through  holes  2,  4  and  5  and 
the  z„  axis  is  in  the  plane  passing  through  holes  1,  3  and  5.  The  positive  directions  of 
the  axes  are  as  shown  and  the  axes  form  a  right-handed  orthogonal  system.  The  axes 
convention  is  the  same  as  that  used  by  Fairlie  (Reference  9).  The  tunnel  coordinate 
system,  xT  yTzT,  introduced  in  Section  2.1,  is  also  shown  in  Figure  14.  The  relationship 
between  the  two  coordinate  systems  is  given  in  Section  3.3. 

The  probe  can  be  used  in  either  the  stationary  or  the  null-sensing  mode.  In  the 
stationary  mode,  the  probe  is  given  a  known  orientation  with  respect  to  the  wind  tunnel 
and  by  measuring  the  pressures  in  the  six  holes  the  velocity  VY  and  the  flow  angles  0Y 
and  <(>Y  can  be  determined  using  calibration  relationships.  Velocities  and  flow  angle's 
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Figure  14.  Six-hole  probe  showing  velocity  vector,  flow  angles  and  coordinate  systems 
(based  on  a  diagram  given  in  Reference  8). 


relative  to  the  tunnel  coordinate  system  can  then  be  determined  (Section  3.5).  In  the 
null-sensing  mode,  the  probe  is  oriented  so  that  the  pressures  in  holes  1  to  4  are  equal, 
which  signifies  that  the  velocity  vector  is  parallel  to  the  probe  longitudinal  axis.  The 
magnitude  of  the  velocity  vector  can  then  be  determined  in  the  conventional  manner 
using  the  total  and  static  pressures  measured  in  holes  5  and  6  respectively,  i.e. 
Vy  =  [2(p5  -  p6)/p]05,  and  the  direction  of  the  velocity  vector  can  be  determined  by 
measuring  the  orientation  of  the  probe  relative  to  the  tunnel  coordinate  system.  The 
null-sensing  mode  of  operation  is  cumbersome  and  time  consuming,  whereas  the 
stationary  mode  is  straightforward  and  it  was  the  method  used  for  the  current 
measurements. 

Melbourne  &  Griss  (Reference  8)  calibrated  the  probe  in  the  LSWT  for  free-stream 
velocities  of  100  and  200  ft/s  (30.5  and  61.0  m/s  respectively)  and  obtained  expressions 
relating  Vy,  0Y  and  cf)Y  to  measured  pressures,  p,  to  p6.  From  their  calibration  measure¬ 
ments,  they  found  that 
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0Y=  1.983  R3- 2.105  R2  +  16.917  R 


for  a  free-stream  velocity  of  100  ft/s  (30.5  m/s)  and 


0Y  =  -  0.050  R3  +  2.686  R2  +  14.013  R 


(7) 


(8) 


for  a  free-stream  velocity  of  200  ft/s  (61 .0  m/s).  In  these  two  expressions,  R  is  given  by 

.  RH(Cp,-cp,)Mcp.,-cJT  (9) 

where  the  pressure  coefficients,  (J  ,  (n  =  1  to  4)  are  given  by 


C  =  -  Pn  Prcf  - 
Pn  0.5  p  V2 


(10) 


Melbourne  &  Griss  also  found  that 

C  —  C  =  (  ~  Pa ) 

Ps  *  0.5  p  V2 


=  -  0.000493  0Y  +  0.004472  0Y  +  1 .003  (11) 


for  a  free-stream  velocity  of  100  ft/s  (30.5  m/s)  and 

Cp5  -  CP(  =  -QP5  =  -  0.000535  e;  +  0.005854  0Y  +  1 .014 


(12) 


for  a  free-stream  velocity  of  200  ft/s  (61.0  m/s).  Equations  7,  8,  11  and  12  were 
determined  using  the  method  of  least  squares  to  fit  data. 

The  flow  angle,  <J)Y  ,  is  related  to  pressure  coefficients  by  the  following  relationship: 


(j)Y  =  tan  1 


(13) 


If  0Y  is  very  small,  the  values  of  Cp  ,  Cp  ,  Cp  and  Cp  will  be  approximately  the  same  so 
that  equation  13  will  be  ill-conditioned  and  may  break  down. 
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In  the  above  calibration  relationships  (equations  7  to  13),  VY,  0Y  and  (])Y  are  not  given 
explicitly  in  terms  of  measured  pressures,  p,  to  p6,  and  it  was  necessary  to  determine  the 
three  variables  iteratively.  The  procedure  used  is  discussed  in  Section  3.5. 

In  the  current  investigation,  the  probe  calibration  applicable  to  a  free-stream  velocity  of 
100  ft/s  (30.5  m/s)  was  only  used  for  flow-angle  measurements  at  UN0M  =  30  m/s  and 
similarly  the  calibration  applicable  to  a  velocity  of  200  ft/s  (61.0  m/s)  was  only  used  for 
measurements  at  UN0M  =  60  m/s. 

3.3  Transformation  of  Coordinates 

The  orientation  of  one  set  of  axes  relative  to  any  other  set  can  be  given  by  three  angles, 
which  are  the  consecutive  rotations  about  the  (right-handed)  axes  z,  y  and  x  in  that  order 
that  carry  one  set  of  axes  into  coincidence  with  the  other  (Reference  9).  In  Figure  14, 
the  angles  that  rotate  the  tunnel  axes,  xT,  yT  and  zT,  into  coincidence  with  the  body  axes, 
xb,  yB  and  zB,  are  referred  to  as  the  probe  yaw  angle,  \|/,  the  pitch  angle,  0,  and  the  roll 
angle,  (j).  The  positive  directions  of  \\f,  0  and  (j)  conform  to  the  conventional  right-hand 
screw  rule  about  the  relevant  axes,  as  shown  in  Figure  14.  For  the  orientation  of  the 
tunnel  and  body  axes  shown  in  Figure  14,  \|/  is  180.0°,  0  is  0.0°  and  <})  is  0.0°. 

The  velocity  vector,  VY,  shown  in  Figure  14,  can  beresolved  into  three  components  of 
velocity  in  the  body  coordinate  system  to  give  UB,  VB  and  WB,  corresponding  to  the  xB, 
yB  and  zB  directions  respectively.  For  arbitrary j/alues  o£\|/,  0  and  <|>,  components  of 
velocity  in  the  tunnel  coordinate  system,  UT,  VT  and^Wj.,  can  be  determined  from 
components  of  velocity  in  the  body  coordinate  system,  UB,  VB  and  WB,  using 

"  UT  1  [  UB  ' 

v,  =  M„  V,  (14) 

.  WT  J  W,  . 

where  I  X,L  is  a  matrix  of  direction  cosines  that  is  used  to  transform  velocities  from  the 

L  J  IB  J  I 

body  coordinate  system  to  the  tunnel  coordinate  system.  The  matrix  LA]  TB  is  given  by 

sin  (j)  sin  0  cos  \\f  cos  $  sin  0  cos 
cos  0  cos  \|/  ,  .  .  x  . 

-  cos  (j)  sin  \|/  +  sin  (p  sin  \|/ 

sin  <b  sin  0  sin  w  cos  <b  sin  0  sin  \(/  -Y 

cos  0  sin  \|/  (15> 

+  cos  cf)  cos  \| /  -  sin  (p  cos  \j/ 

-  sin  0  sin  ((>  cos  0  cos  ((>  cos  0 

which  is  the  transpose  of  |Aj  bt,  given  for  example  in  Reference  9. 
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3.4  Traversing  Mechanism 

A  specially  designed  traversing  mechanism  was  used  to  support  the  six-hole  probe  and 
to  position  it  at  different  vertical  locations  in  the  tunnel  when  measuring  velocities  and 
flow  angles.  A  photograph  of  the  mechanism  mounted  in  the  tunnel  is  given  in 
Figure  15  and  a  plan  cross  section  assembly  drawing  of  the  mechanism  showing  the 
main  features  is  given  in  Figure  16.  Photographs  of  components  of  the  mechanism  are 
given  in  Figure  17.  The  probe  is  mounted  on  a  carriage  which  can  move  up  or  down 
(into  and  out  of  the  paper  in  Figure  16)  between  two  linear-motion  rails,  and  the 
carriage  is  connected  to  the  linear-motion  rails  using  three  low-friction  linear-motion 
bearings  -see  Figure  17  (b).  A  stepping  motor  is  connected  to  a  lead  screw  which 
passes  through  the  carriage  and  instructions  can  be  sent  to  the  motor  from  a  computer  to 
rotate  the  lead  screw  and  move  the  probe  to  a  specified  zT  location  -see  Figure  17  (c). 
Limit  switches  are  installed  near  the  top  and  near  the  bottom  of  the  downstream  linear- 
motion  rail  to  prevent  the  carriage  moving  beyond  acceptable  lower  and  upper  limits  of 
travel  -see  Figure  17.  An  electric  brake,  which  is  activated  when  the  switch  to  the 
stepping  motor  is  off  or  when  the  power  to  the  motor  fails,  was  installed  on  the  lead 
screw  to  prevent  it  from  rotating  under  these  conditions  -see  Figure  17  (a).  Without  the 
brake,  the  weight  of  the  carriage  and  the  components  attached  would  cause  the  lead 
screw  to  rotate  and  the  carriage  would  move  down  when  there  was  no  power  to  the 
motor.  A  dummy  probe  is  attached  to  the  traversing  mechanism  to  counteract  possible 
twisting  of  the  mechanism  that  would  arise  from  aerodynamic  loads  on  a  single  probe. 
The  active  and  dummy  probes  can  be  swapped  as  required  by  simply  transferring  the 
probes  from  one  collet  to  the  other  (Figure  16),  without  changing  the  yaw,  pitch  and  roll 
angles  of  the  probes.  The  horizontal  arm  holding  the  probes  is  interchangeable  with 
arms  of  different  lengths,  thereby  enabling  measurements  to  be  taken  at  different  yT 
locations  without  changing  the  lateral  position  of  the  rig. 

The  traversing  mechanism  was  bolted  to  two  rectangular  plates,  one  of  which  was 
bolted  to  the  floor  and  the  other  to  the  roof  of  the  test  section.  The  plates  were  0.4  m 
long  and  1.6  m  wide  and  occupied  most  of  the  width  of  the  floor  (or  roof)  not  counting 
the  fillets  (the  lower  plate  is  shown  in  Figure  15).  The  lower  plate  was  installed  first 
and  it  was  fixed  to  the  floor  of  the  test  section  so  that  its  trailing  edge  was  parallel  to  the 
yT  direction  and  so  that  the  tip  of  the  six-hole  probe  would  be  located  at  the  required  xT 
position  when  the  mechanism  was  bolted  between  the  plates.  The  lower  plate  was 
aligned  accurately  by  first  ensuring  that  the  midpoint  of  its  trailing  edge  was  positioned 
on  the  longitudinal  centreline  of  the  floor  of  the  tunnel.  If  necessary,  the  yaw  angle  of 
the  plate  was  then  adjusted  slightly  so  that  the  distances  from  each  of  the  two 
extremities  of  the  trailing  edge  to  an  upstream  or  downstream  reference  mark  (at  least 
2.0  m  from  the  trailing  edge)  on  the  longitudinal  centreline  of  the  floor  of  the  tunnel 
were  the  same.  Using  this  technique,  the  error  in  the  setting  of  the  yaw  angle  of  the 
plate  could  be  kept  to  less  than  about  ±0.1°.  The  alignment  of  the  trailing  edge  of  the 
plate  with  the  yT  direction  was  critical  since  probe  yaw  angles  were  referenced  to  the 
plate  trailing  edge.  The  fore/aft  location  of  the  lower  plate  did  not  have  to  be  precise, 
since  the  probe  could  be  adjusted  axially  (by  up  to  about  45  cm)  through  the  probe 
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Figure  16.  Assembly  drawing  of  the  traversing  mechanism  and  probes. 


(b)  carriage 


(c)  bottom  part  showing  stepping  motor 


Figure  17.  Components  of  the  traversing  mechanism. 


(a)  top  part  showing  electric  brake 
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collet  (Figure  16).  A  special  alignment  frame,  shown  in  Figure  18,  was  used  to  ensure 
that  the  upper  plate  was  located  directly  above  the  lower  plate.  The  frame  was  bolted 
firmly  to  the  lower  plate,  using  dowel  pins  for  alignment.  The  top  plate  was  then 
mounted  on  four  blocks  attached  to  the  frame,  using  dowel  pins  for  alignment.  Four  set 
screws  on  the  frame  were  used  to  lift  the  blocks  and  the  top  plate  until  the  plate  came 
into  contact  with  the  roof  of  the  test  section.  Guide  pins  minimized  lateral  movement  of 
the  top  plate  while  it  was  being  raised.  The  top  plate  was  then  bolted  to  the  roof  of  the 
test  section  and  the  frame  was  removed  from  the  tunnel.  The  traversing  mechanism  was 
then  positioned  on  the  lower  plate,  using  a  guide  bar  at  the  foot  of  the  mechanism 
(Figure  16)  to  ensure  that  the  base  of  the  mechanism  was  not  twisted  relative  to  the 
trailing  edge  of  the  lower  plate.  The  traversing  mechanism  was  then  bolted  to  the  lower 
plate  so  that  the  probe  tip  was  located  at  the  required  yT  position.  The  two  vertical 
linear-motion  rails  on  the  mechanism  (Figure  16)  were  positioned  so  that  they  were  less 
than  0.05°  from  the  vertical,  as  measured  with  an  inclinometer  on  two  adjacent  sides  of 
each  rail,  and  the  mechanism  was  then  bolted  to  the  upper  plate. 

There  is  provision  on  the  mechanism  to  change  the  orientation  of  the  six-hole  probe 
(and  the  dummy  probe).  The  yaw  angle,  \| /,  pitch  angle,  0,  and  roll  angle,  <>,  can  all  be 
altered  independently  (see  Figure  16  for  adjustment  mechanisms).  On  the  mechanism, 
V  could  be  varied  by  about  ±15°,  6  by  360°  and  <(>  by  360°.  It  would  have  been 
relatively  difficult  to  set  values  of  \)/  to  the  accuracy  needed  when  the  probe  longitudinal 
axis  was  horizontal  (i.e.  0  =  0°),  but  by  rotating  the  pitch  shaft  (Figure  16)  so  that  the 
probe  longitudinal  axis  was  temporarily  approximately  vertical,  and  then  adjusting  the 
yaw  screws  to  set  the  probe  vertical  in  the  yT-zT  plane,  using  an  inclinometer  on  the 
probe  body,  and  then  returning  the  probe  to  the  horizontal  position,  values  of  \ /  could  be 
set  to  0.0°  to  an  accuracy  of  about  ±0.02°.  An  inclinometer  was  also  used  to  set  values 
of  0  to  an  accuracy  of  about  ±0.02°.  Values  of  <f»  were  set  by  rolling  the  probe  until  the 
scribed  reference  line  on  the  body  of  the  probe  (Figure  14)  matched  the  appropriate 
scribed  line  on  the  graduated  scale  on  the  collet  (Figure  13).  A  microscope  was  used 
when  matching  the  scribed  lines  and  values  of  <|)  were  set  to  an  accuracy  of  about  ±0.2°. 

The  traversing  mechanism  is  a  precision  device  and  was  manufactured  to  high  accuracy. 
However,  as  with  any  practical  device,  the  individual  components  on  the  mechanism  all 
had  manufacturing  tolerances.  Checks  were  done  and  it  was  found  that  the  possible 
error  in  the  setting  of  the  yaw  angle  of  the  six-hole  probe  due  to  these  effects  could  be 
up  to  ±0.1  °.  This  would  occur  if  components  had  been  manufactured  to  the  limits  of 
their  tolerances  (either  maximum  or  minimum  limits)  and  then  assembled  so  that  the 
effects  of  the  inbuilt  imperfections  were  maximized,  i.e.  the  worst  possible  case. 

It  was  found  that  the  engraved  scales  on  both  of  the  probe  collets  (Figure  13)  had  been 
rolled  slightly  from  the  true  zero  position  and  these  offsets  have  to  be  allowed  for  when 
setting  probe  roll  angles.  The  engraved  scale  on  the  port  collet  has  a  zero  error  of 
+2.10°  and  the  scale  on  the  starboard  collet  has  a  zero  error  of  +1.35°,  i.e.  the  engraved 
scales  are  both  rolled  clockwise  when  looking  in  the  direction  from  the  rear  of  the  probe 
to  the  hemispherical  tip  of  the  probe  (Figure  13). 
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Figure  1 8.  Alignment  frame  used  when  locating  the  upper  plate  of  the  traversing 
mechanism  on  the  test  section  roof  with  respect  to  the  lower  plate  on  the  floor. 


The  movement  of  the  probe  was  controlled  using  program  WINDRIG  located  on  a  PC. 
Initially  the  user  must  enter  the  current  zT  position  of  the  probe  relative  to  the  floor  of 
the  test  section.  This  initial  setting  is  used  as  a  reference  for  subsequent  positioning  of 
the  probe.  The  user  must  then  enter  the  lower  and  upper  zT  limits  corresponding  to  the 
permissible  probe  movement.  These  limits  are  independent  of  the  limits  set  by  the  limit 
switches  (see  above)  and  would  be  used,  for  example,  if  the  probe  was  being  traversed 
between  the  comer  fillets  in  the  octagonal  test  section.  To  minimize  the  possibility  of 
damaging  the  probe  and  the  traversing  mechanism,  the  program  has  been  written  so  that 
it  is  not  possible  to  move  the  probe  using  the  program  unless  these  limits  have  been 


39 


DSTO-TR- 1 384 


entered.  Finally,  the  user  must  enter  the  destination  value  of  z,  and  the  mechanism  then 
moves  the  probe  as  required.  The  smallest  permissible  probe  movement  from  a  given 
location  is  0.1  mm  and  the  probe  can  be  positioned  to  an  accuracy  of  0.1  mm  relative  to 
its  initial  setting. 


3.5  Procedure  Used  to  Measure  Flow  Angles 


When  measuring  velocities  and  flow  angles,  the  probe  was  positioned  in  the  tunnel  so 
that  the  origin  of  the  body  coordinate  system,  x„  y„  zB  (Figure  14),  was  located  at  the 
grid  point  at  which  measurements  were  to  be  taken.  For  all  of  the  current 
measurements,  the  probe  orientation  angles,  \| /,  0  and  {>,  were  set  at  180.0°,  -1.0°  and 
0.0°  respectively  (Figure  14).  The  probe  vertical  movement  was  restricted  by  the  brake 
(at  the  top  of  the  mechanism)  and  by  the  stepping  motor  (at  the  bottom  of  the 
mechanism)  and  it  was  necessary  to  set  9  at  -1.0°  to  enable  measurements  to  be  taken  at 
both  the  upper  and  lower  grid  points  at  given  yT  locations. 

To  determine  the  direction  of  the  flow  (Figure  14),  it  was  first  necessary  to  measure 
probe  pressures,  p,  to  p6.  These  pressures  were  measured  using  the  same  electronic 
manometer  system  and  data  sampling  format  that  were  used  to  measure  the  pressures 
for  the  mean  velocities,  as  described  in  Section  2.4.  For  the  flow  angle  measurements,  a 
PSI  Electronic  Pressure  Scanner  having  a  full-scale  range  of  ±2490  Pa  (±10.0  inches  of 
water)  was  used. 

Once  p,  to  p6  had  been  measured,  VY,  0Y  and  <J)Y  were  determined  using  the  probe 
calibration  (equations  7  to  13).  As  indicated  in  Section  3.2,  the  calibration  relationships 
do  not  give  VY,  0Y  and  <j>Y  explicitly  in  terms  of  measured  pressures,  p,  to  p6,  and  it  was 
necessary  to  determine  the  three  variables  iteratively.  The  following  procedure  was 
used. 

1 .  Asa  first  approximation,  the  value  of  VY  was  set  below  its  known  value  (e.g.  VY  was 
initially  set  at  27  m/s  when  the  free-stream  velocity  at  which  the  measurements  were 
taken  was  known  to  be  about  30  m/s). 

2.  Using  this  first  approximation  of  VY,  first  approximations  of  C  (n  =  1  to  4),  R  and 
0Y  were  made  using  equations  10,  9  and  7  (or  8)  respectively.  Pn 

3.  Using  the  first  approximation  of  0Y,  a  second  approximation  of  VY  was  made  using 
equation  1 1  (or  12). 

4.  Successive  approximations  of  VY  were  made  and  steps  1  to  3  were  repeated  until  the 

velocities  converged  to  within  0.001  m/s.  At  this  stage,  Vy,  0Y,  C  C  ,  C  and  C 
were  determined  to  an  acceptable  accuracy.  Pl  Pl  Pj  P4 

5.  Once  values  of  Cp  were  known,  <))Y  was  calculated  using  equation  13. 
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Figure  19.  Velocities  and  flow  angles  in  the  tunnel  coordinate  system. 

VY  was  resolved  into  UB,  VB  and  WB,  the  components  of  velocity  in  the  xH,  yB  and  zB 
directions  respectively  in  the  body  coordinate  system,  and  these  velocities  were 
transposed,  using  equation  14,  into  UT,  VT  and  WT,  the  components  of  velocity  in  the  xT, 
yT  and  zT  directions  respectively  in  the  tunnel  coordinate  system.  Horizontal  and 
vertical  flow  angles,  £H  and  respectively,  relative  to  the  xT  axis  as  shown  in  Figure  19, 
were  then  calculated.  All  iterative  calculations  and  data  reduction  were  carried  out 
using  the  program  ANGLE  located  on  a  PC. 

3.6  Corrections  to  Flow  Angles  to  Allow  for  Probe  Support  Interference 

The  traversing  mechanism  used  to  position  the  six-hole  probe  in  the  tunnel  when 
measuring  flow  angles  was  made  to  be  reasonably  rigid  (Figures  15  and  16).  The 
vertical  column  on  the  mechanism  is  46  mm  wide  and  extends  from  the  floor  to  the  roof 
of  the  test  section,  so  that  the  ratio  of  the  cross-sectional  area  of  the  column  to  the  total 
transverse  cross-sectional  area  of  the  test  section  is  0.019:1,  i.e.  1.9%.  Although  such  a 
value  of  blockage  is  significantly  less  than  the  value  7.5%,  quoted  in  Reference  10  as 
being  the  upper  limit  of  blockage  acceptable  for  wind-tunnel  testing,  interference 
caused  by  the  column  nevertheless  had  a  noticeable  effect  on  measured  flow  angles. 
When  the  six-hole  probe  was  supported  on  the  port  side  of  the  vertical  column,  flow 
angles  measured  at  given  zT  locations  on  the  vertical  centreplane  in  the  test  section  were 
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found  to  be  noticeably  different  (see  below)  from  those  measured  at  the  same  locations 
when  the  probe  was  supported  on  the  starboard  side  of  the  vertical  column.  Note  that 
the  vertical  column  on  the  traversing  mechanism  had  to  be  moved  from  one  side  of  the 
test  section  to  the  other  side  to  obtain  the  two  sets  of  measurements. 

Experiments  were  undertaken  to  determine  the  corrections  to  be  applied  to  measured 
flow  angles  to  allow  for  the  effects  of  the  interference  caused  by  the  traversing 
mechanism.  To  determine  a  typical  interference  correction,  the  six-hole  probe  was 
mounted  in  a  specially-designed  thin  supporting  column  (width  of  thin  column  is 
16  mm),  so  that  the  tip  of  the  probe  was  directly  upstream  of  the  column  and  the  probe 
yaw,  pitch  and  roll  angles  were_set  close  to  their  values  for  normal  operation.  Flow 
angles  were  then  measured  for  UNOM  =  30  and  60  m/s.  The  traversing  mechanism  was 
not  in  the  test  section  for  these  measurements.  The  traversing  mechanism  was  then 
positioned  in  the  test  section  on  either  the  port  side  of  the  probe  or  on  the  starboard  side 
of  the  probe,  but  not  actually  attached  to  the  probe,  and  flow  angles  were  measured 
again.  The  vertical  column  on  the  traversing  mechanism  was  located  at  the  same  lateral 
distance  from  the  probe  as  it  would  be  during  normal  operation,  when  the  probe  was 
attached  to  the  mechanism.  A  typical  experimental  setup  is  shown  in  Figure  20  for  the 
case  when  the  thin  column  and  the  traversing  mechanism  are  both  in  the  test  section. 
The  interference  corrections  were  then  determined  as  the  differences  in  corresponding 
flow  angles  measured  both  with  the  traversing  mechanism  in  the  tunnel  and  with  the 
mechanism  out  of  the  tunnel. 


When  the  six-hole  probe  was  mounted  on  the  thin  column  to  assess  probe  support 
interference,  it  did  not  have  to  be  set  precisely  at  the  yaw,  pitch  and  roll  angles  used 
when  the  probe  was  attached  to  the  traversing  mechanism  during  normal  operation. 
When  determining  interference  corrections,  the  important  measurements  were  the 
differences  in  flow  angles  measured  when  the  traversing  mechanism  was  in  the  tunnel 
and  when  it  was  not  in  the  tunnel,  rather  than  the  absolute  values  of  flow  angles. 

Interference  corrections  were  determined  in  the  xT  =  0.0  m  transverse  plane  of  the  test 
section  for  y^W  =  0.0,  ±0.15  and  ±0.3,  and  z/H  =  0.0  and  ±0.3,  where  W  =  2743  mm 
and  H  =  2134  mm  are  the  width  and  height  respectively  of  the  test  section  of  the  LSWT 
(Figure  2).  It  was  necessary  to  use  three  different  thin  supporting  columns  to  support 
the  six-hole  probe  at  the  three  different  zT  locations  in  the  test  section.  The  three 
columns  were  essentially  the  same,  but  had  the  probe  attachment  sleeve  located  at  a 
different  height  (zT  location).  Altogether,  interference  corrections  were  obtained  for  15 
locations,  and  corrections  at  other  locations  at  which  flow  angles  were  measured  were 
obtained  by  linear  interpolation._The  interference  corrections  varied  depending  on  the 
location  in  the  test  section.  For  UNOM  =  30  and  60  m/s,  average  values  of  corrections  for 
were  about  0.4  and  those  for  were  about  0.0°.  The  corrections  determined  for  the 
xT  =  0.0  m  transverse  plane  were  applied  to  flow  angles  measured  in  this  plane  as  well 
as  in  the  xT  =  -2.0,  -1.0  and  1.0  m  transverse  planes.  The  overlapping  results  at  the 
vertical  centreplane  were  effectively  the  same  after  the  corrections  had  been  applied. 
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Traversing  mechanism 
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to  traversing 
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Figure  20.  Traversing  mechanism  and  thin  column  (supporting  the  six-hole  probe) 
mounted  in  the  test  section  of  the  low-speed  wind  tunnel. 
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3.7  Analysis  of  Flow  Angles 

Contours  of  for  UNOM  =  30  m/s  for  xT  =  -2.0,  -1.0,  0.0  and  +1.0  m  are  shown  in 
Figure  21.  Corresponding  plots  for  ^  are  shown  in  Figure  22.  Figures  23  and  24  show 
contours  for  UN0M  =  60  m/s.  As  for  the  mean  velocities,  Figures  21  to  24  each  indicate 
how  the  flow  angles  change  as  the  flow  moves  downstream.  For  the  two  values  of  UNOM, 
the  general  pattern  of  the  £„  contours  at  the  four  xT  locations  is  similar  (Figures  21  and 
23),  with  little  evidence  of  the  flow  becoming  more  uniform  as  it  moves  downstream, 
and  likewise  for  the  corresponding  ^  contours  (Figures  22  and  24).  In  the  central  50% 
of  the  cross-sectional  area  of  the  test  section,  between  xT  =  -2.0  and  +1.0  m,  for 
Unom  =  30  m/s,  horizontal  flow  angles  varied  from  -0.8°  to  +0.4°  (Figure  21)  and 
vertical  flow  angles  varied  from  -1.0°  to  +0.4°  (Figure  22).  For  UN0M  =  60  m/s, 
horizontal  flow  angles  varied  from  -0.6°  to  +0.8°  (Figure  23)  and  vertical  flow  angles 
varied  from -1.2°  to  +0.4°  (Figure  24).  However,  between  xT  =  -1.0  and  +1.0  m,  for 
Unom  =  30  m/s,  horizontal  flow  angles  varied  from  -0.8°  to  +0.4°  and  vertical  flow 
angles  varied  from  -0.8°  to  +0.4°.  For  UNOM  =  60  m/s,  horizontal  flow  angles  varied 
from  -0.4°  to  +0.8°  and  vertical  flow  angles  varied  from  -1.0°  to  +0.4°.  There  is  a 
pocket  of  flow  in  the  lower  starboard  corner  of  the  test  section  which  has  significantly 
larger  values_of  vertical  flow  angles  which  range  up  to  -1.4°  (Figure  24).  Figures  21 
and  22,  for  UNOM  =  30  m/s,  and  Figures  23  and  24,  for  UN0M  =  60  m/s,  indicate  that 
regions  of  high  or  low  horizontal  flow  angles  are  not  necessarily  associated  with 
corresponding  regions  of  high  or  low  vertical  flow  angles. 

Selected  data  associated  with  the  contours  shown  in  Figures  21  to  24  have  been 
replotted  in  Figure  25  to  indicate  how  horizontal  and  vertical  flow  angles  in  the 
horizontal  and  vertical  centreplanes  vary  in  the  longitudinal  or  xT  direction.  The  method 
of  plotting  the  data  is  similar  to  that  used  in  Section  2.5  for  the  mean  velocities.  Data 
corresponding  to  the  xT-zT  plane  that  passes  through  yT  =  0.0  m  (i.e.  the  vertical  plane 
that  is  aligned  with  and  passes  through  the  tunnel  longitudinal  axis)  are  shown  for 
z/H  =  0.0,  ±0.2  and  ±0.4.  Similarly,  data  corresponding  to  the  xT-yT  plane  that  passes 
through  zT  =  0.0  m  (i.e.  the  horizontal  plane  that  is  aligned  with  and  passes  through  the 
tunnel  longitudinal  axis)  are  shown  for  y/W  =  0.0,  ±0.15,  ±0.30  and  ±0.45.  In 
Figure  25,  data  corresponding  to  the  longitudinal  centreline  of  the  tunnel  are  represented 
by  thick  lines.  Horizontal  and  vertical  flow  angles  for  the  chosen  locations  remain 
approximately  constant  over  the  xT  range. 

3.8  Errors  in  Measured  Flow  Angles 

There  are  small  errors  in  the  measured  horizontal  and  vertical  flow  angles,  and 
respectively,  due  to  (1)  instrumentation  errors  and  (2)  probe  setup  errors.  The  probe 
setup  errors  occur  because  the  six-hole  probe  is  not  set  at  precisely  the  required  yaw, 
pitch  and  roll  angles,  relative  to  the  tunnel  coordinate  system,  when  measuring  flow 
angles. 


44 


DSTO-TR-1384 


3.8.1  Instrumentation  Errors 

Considering  instrumentation  errors,  it  was  necessary  to  measure  air  temperatures  and  air 
pressures  when  determining  flow  angles.  The  sensing  devices  used  to  measure 
temperatures  were  only  accurate  to  within  ±0.5°  C  and  the  indicated  readings  of  the 
scanners  used  to  measure  pressures  were  only  accurate  to  within  +0.05%  of  their  full- 
scale  reading.  Two  different  types  of  scanners  were  used  when  measuring  pressures 
and  their  ranges  were  159  kPa  (23.0  psi)  (used  to  measure  reference  velocities)  and 
±2490  Pa  (±10.0  inches  of  water).  The  possible  errors  in  measured  values  of  £H  and  £v 
corresponding  to  these  instrumentation  errors  were  determined.  For  UN0M  =  30  m/s, 
errors  in  £„  were  about  ±0.08°  and  errors  in  £y  were  about  ±0.08°.  Corresponding 
errors  for  UNOM  =  60  m/s  were  about  ±0.02°  for  £(1  and  about  ±0.02°  for  £v. 

3.8.2  Probe  Setup  Errors 

Considering  probe  setup  errors,  it  is  indicated  in  Section  3.4  that  the  possible  errors  in 
setting  the  yaw,  pitch  and  roll  angles  of  the  six-hole  probe  are  ±0.02°,  ±0.02°  and  ±0.2° 
respectively.  The  ±0.02°  probe  setup  error  given  for  the  yaw  angle  corresponds  to  the 
possible  misalignment  of  the  probe  relative  to  the  traversing  mechanism,  and  not  the 
misalignment  of  the  probe  relative  to  the  tunnel.  It  is  also  indicated  in  Section  3.4  that 
an  additional  error  of  ±0.1°  in  probe  yaw  angle  may  exist  since  the  traversing 
mechanism  itself  can  only  be  set  to  an  accuracy  of  0.1  °  in  yaw  angle,  due  to  difficulties 
in  accurately  positioning  the  lower  plate  which  is  used  to  align  the  mechanism.  It  is 
further  indicated  that  an  additional  error  of  ±0.1°  in  probe  yaw  angle  may  exist  due  to 
possible  twisting  of  the  mechanism  relative  to  the  lower  plate,  resulting  from  imperfect 
manufacture  and  assembly  of  the  mechanism.  Possible  errors  in  measured  values  of  £„ 
and  Ey  corresponding  to  the  probe  setup  errors  given  above  were  determined.  It  was 
found  that  errors  in  £,,  were  about  ±0.02°  and  errors  in  £v  were  about  ±0.05°.  The  setup 
error  of  ±0.1°  in  probe  yaw  angle,  due  to  the  positioning  the  lower  plate,  results  in  an 
0.1°  increase  or  decrease  in  the  possible  error  in  £„,  and  similarly  for  the  setup  error  of 
0.1°  in  probe  yaw  angle  due  to  the  twisting  of  the  mechanism.  These  latter  two  setup 
errors  in  yaw  angle  have  negligible  effect  on  errors  in  measured  values  of  £v. 

3.8.3  Accumulation  of  Errors 

All  possible  errors  in  £„  and  ^  given  above  are  additive,  i.e.  the  errors  in  £„  and  ^  can 
be  as  large  as  ±0.22°  and  ±0.05°  respectively.  However,  since  the  lower  plate  was  not 
moved  during  any  set  of  measurements  at  a  given  xT  location,  any  errors  due  to  the 
positioning  of  the  lower  plate  will  affect  all  measured  values  of  £„  equally  at  that  xT 
location,  i.  e.  all  measured  values  of  £„  for  a  particular  cross  section  of  the  flow  may  be 
high  or  low  by  up  to  0.1°  due  to  this  effect.  Similarly,  since  the  traversing  mechanism 
was  not  dismantled  after  initial  assembly,  any  errors  due  to  the  twisting  of  the 
mechanism  will  affect  all  measured  values  of  £„  equally  for  the  entire  set  of  angle 
measurements,  i.  e.  all  measured  values  of  £,,  for  the  entire  experimental  program  may 
be  high  or  low  by  up  to  an  additional  0.1  °. 


45 


£h  deg 

Unom  =  60  m/s 
xT  =  -2.0  m 


£h  deg 

UNom  =  60  m/s 
xT  =  -1.0  m 


1.0,  0.0  and  1.0  m. 


51 


I 


mtM+  twin 


2 


(c)  Horizontal  flow  angles,  UN0M  =  60  m/s. 


xT-zT  plane,  yT  =  0.0  m 


£v 

(deg) 


A 

w 

— 0.45 
— O-  0.30 
— n-  o.i5 

— O — 0.15 
— B—  -0.30 
-V-  -0.45 
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4.  Turbulence-Intensity  Measurements 

4.1  Test  Schedule 

Turbulence  intensities,  (^)°7Uon,  (^)°7U00  and  (^)°'7U0n,  were  measured  in  the  test 
section  at  the  31  grid  points  shown  in  Figure  26  (a)  for  xT  =  0.0  m  as  well  as  at  the  15 
grid  points  shown  in  Figure  26  (b)  for  xT  =  -2.0,  -1.0  and  +1.0  m.  Intensities  were 
measured  for  UNOM  =  30  and  60  m/s.  The  symbols  u,  v  and  w  denote  the  fluctuating 
components  of  velocity  in  the  xT,  yT  and  z^  directions  respectively.  (Note  that  in  a 
fluctuating  flow,  U,  the  instantaneous  velocity  in  the  xT  direction,  is  related  to  the  mean 
and  fluctuating  components  of  velocity  in  that  direction  by  U  =U  +  u,  and  likewise  for 
the  yT  and  zT  directions). 


4.2  Modified  Cross-Tunnel  Rake 


For  the  intensity  measurements,  the  cross-tunnel  rake  (Section  2.3)  was  modified  by 
removing  some  of  the  pressure  probes  and  replacing  them  with  posts,  which  were  used 
to  attach  a  hot-wire  probe. 

43  Procedure  Used  to  Measure  Turbulence  Intensities 

The  turbulence  intensities  were  measured  using  the  PSL  constant-temperature  hot-wire 
anemometer  system  described  in  Reference  1 1.  The  system  includes  four  anemometers, 
a  matching  box,  a  dynamic  calibrator,  four  low-pass  filters,  a  data-acquisition  system 
(based  on  a  PC)  and  software  for  calibrating  the  hot  wires  and  for  acquiring  and 
reducing  the  data.  For  the  current  measurements,  a  DANTEC-  55P51  crossed^ire 
probe  and  corresponding  leads  were  used,  but  the  platinum-plated  tungsten  wires  on  the 
probe  were  removed  and  replaced  with  Wollaston  wires.  Although  platinum-plated 
tungsten  wires  are  stronger,  Wollaston  wires  were  used  because  the  existing  hot-wire 
anemometers  had  been  optimized  for  use  with  this  type  of  wire.  Wollaston  wires  also 
have  a  more  stable  temperature  calibration  than  platinum-plated  tungsten  wires  (see 
Reference  11). 

The  procedures  used  to  replace  the  hot  wires,  set  the  operating  resistances  of  the  wires, 
set  the  frequency  responses  of  the  anemometers,  electronically  match  the  wires, 
calibrate  the  wires  using  the  dynamic  calibrator,  sample  the  wire  voltages  using  the 
data-acquisition  system  and  reduce  the  data  to  obtain  the  required  turbulence  terms,  are 
described  in  Reference  1 1. 

On  the  modified  cross-tunnel  rake  (Section  4.2),  only  one  hot-wire  probe  was  used  for 
the  intensity  measurements  and  it  was  necessary  to  roll  the  probe  on  its  post  or  move  the 
probe  from  post  to  post  to  enable  the  complete  set  of  measurements  to  be  taken.  After 

2  DANTEC  Elektronic,  Tonsbakken  16-18,  DK-2740,  Skovlunde,  Denmark. 
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Figure  26.  Grid  points  at  which  turbulence  intensities  (unfilled  plus  filled  circles)  and 
spectra  (filled  circles)  were  measured,  (a)  xT  =  0.0  m,  (b)  xT  =  -2.0,  -1 .0  and  1 .0  m. 


DSTO-TR- 1 384 


measuring  intensities  for  the  u  and  v  components  of  the  turbulence,  the  probe  was  rolled 
90°  to  measure  intensities  for  the  u  and  w  components  of  the  turbulence.  At  each  probe 
location,  the  u-component  intensities  were  determined  by  averaging  the  two  measured 
intensities  corresponding  to  the  two  roll  anglesjrf  the  probe.  The  difference  between 
the  results  for  the  two  roll  angles  was  small  for  UN0M  =  30  and  60  m/s,  being  typically  of 
the  order  of  0.02  and  0.06  percentage  points  respectively.  The  tunnel  had  to  be  stopped 
each  time  the  probe  was  rolled  or  relocated. 

When  measuring  intensities,  each  hot-wire  voltage  was  sampled  30  000  times  at  a 
sampling  frequency  of  200  Hz  and  the  low-pass  filters  were  set  at  10  kHz.  These  samp¬ 
ling  parameters  have  been  shown  to  be  satisfactory  for  the  measurement  of  intensities  in 
the  free-stream  of  the  LSWT  (Reference  11).  Since  the  root-mean-square  intensities 
were  known  to  be  much  less  than  10%  at  all  locations  in  the  flow,  hot-wire  non-linearity 
was  neglected  and  intensities  were  evaluated  by  processing  measured  hot-wire  voltages 
rather  than  by  processing  velocities  determined  from  voltages  (Reference  12).  The  hot 
wires  were  calibrated  at  both  UN0M  =  30  and  60  m/s  and  hot-wire  sensitivities  were 
evaluated  at  these  two  velocities.  The  turbulence  intensities  were  then  determined  using 
these  sensitivities.  Since  sensitivities  were  only  required  at  either  UN0M  =  30  or  60  m/s, 
it  was  not  necessary  to  carry  out  a  complete  dynamic  calibration  to  obtain  calibration 
equations  covering  a  range  of  free-stream  velocities. 

For  each  position  of  the  hot-wire  probe,  valuesof  (u2) 0  \  (v2)0'  and  (w1)"5  were  measured 
and  at  the  same  time  corresponding_values  of  URI:|  were  recorded.  Using  these  values  of 
URfj:,  longitudinal  mean  velocities,  U,  at  each  probe  location  were  then  determined  from 
the  mean-velocity  measurements  (Section  2).  In  general,_thejvalues  of  U  varied  slightly 
across  the  cross  section.  Turbulence  intensities,  (u2)"7U,  (v2)°7U  and  (w5)0'/!!,  were 
then  computed  at  each  point  in  the  grid,  using  the  appropriate  value  of  U  for  that 
location  as  the  normalizing  velocity.  The_calculated  intensities  were  then  corrected  very 
slightly  so  that  they  were  defined  using  U()  n,  the  mean  velocity  at  the  centre  of  the  test 
section. 


4.4  Analysis  of  Turbulence  Intensities 

Contours  of  (u  )  /U0()  for  UN0M  =  30  m/s  at  Xj.  =  -2.0^-ri  .0,_0.0  and  +1 .0  m  are  shown  in 
Figure  27.  Corresponding  plots  for  (v2)°7u,,n  and  (w2)°7Uon  are  shown  in  Figures  28 
and  29  respectively.  Figures  30,  31  and  32  show  contours  for  UN0M  =  60  m/s.  Each  of 
Figures  27  to  32  have  been  plotted  to  indicate  how  the  contours  change  as  the  flow 
move^downstream,  i.e.  xT  changes.  For  the  two  values  of  UN0M,  the  contours  of 
(u  )  /Uon  at  the  four  xT  locations  (Figures  27  and  30)  have  the  same  general  features,  as 
far  as  it  is  possible  to  tell  from  the  limited  data  at  xT  =  -2.0,  -1.0  and  +1.0  m,  and 
likewise  for  the  corresponding  contours  of  ( v1)0 5/U() 0  (Figures  28  and  31  respectively) 
and  (w2)  /Uon  (Figures  29  and  32  respectively).  This  indicates  that  regions  of  a  given 
turbulence  intensity  are  convected  downstream  with  little  lateral  spreading  and  there 
does  not  appear  to  be  any  significant  decay  of  the  turbulence  as  the  flow  moves 
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downstream.  The  most  important  result  is  the  value  of  (u2)°7U00,  which  is  a  measure  of 
the  turbulence  level  in  the  tunnel.  In  the  central_50%  of  the  cross-sectional  area  of  the 
test  section,  for  UN0M  =  30  and  60  m/s,  (u2)°7U00  is  generally  below  about  04%,  at 
xT  =  -2.0,  -1.0,  0.0  and  +1.0  m,  whereas  corresponding  values  of  (v2)°'7U00  and 
( w7) 0  5/U0  n  are  generally  below  about  0.7%  (Figures  27  to  32).  The  intensity  levels  for 
the  longitudinal  (u)  component  of  the  turbulence  are  about  two  thirds  those  for  the 
lateral  (v  and  w)  components.  This  non-isotropic  behaviour  is  common  in  wind-tunnel 
flows  since  the  longitudinal  vorticity  component  is  reduced  by  a  greater  amount  than  the 
lateral  components  by  the  contraction  (Reference  13).  The  strongest  gradients  of 
intensity  occur  on  the  starboard  side  of  the  test  section  where  the  intensities  attain  their 
maximum  values. 

4.5  Errors  in  Measured  Turbulence  Intensities 

There  are  errors  in  the  measured  turbulence  intensities  due  to  instrumentation  errors  and 
also  due  to  the  fact  that  the  calibration  of  the  hot  wires  may  have  drifted  slightly 
between  the  start  and  finish  of  a  set  of  measurements.  Instantaneous  velocities 
measured  by  the  hot  wires  were  accurate  to  within  about  ±1.5%  due  to  the  drifting^of 
the  calibration  ofthe  hot  wires.  It  was  found  that,  for  UN0M  =  30  m/s,  errors  in  (u2)°7Uno, 
(v2)°',/U00  and  (w2)  °5/U00  were  about  ±0.011,  ±0.018^  andL±0.0 18  percentage  points 
respectively.  For  example,  if  the  measured  value  of  (u2)°  7U00  at  30  m/s  was  0.5%  then 
its  value  could  range  from  0.51 1%  to  0.489%.  Corresponding  errors  for  UNOM  =  60  m/s 
were  about  ±0.007,  ±0.013  and  ±0.013  percentage  points  respectively. 


5.  Spectral  Measurements 


5.1  Test  Schedule 

Spectra  for  the  u-,  v-  and  w-components  of  the  turbulence  were  measured  in  the  test 
section  at  the  5  grid  points  shown  in  Figure  26  (a)  for  xT  =  0.0  m  as  well  as  at  the  5  grid 
points  shown  in  Figure  26  (b)  for  xT  =  -2.0,  -1.0  and  +1.0  m.  Spectra  were  measured 
for  UN0M  =  30  and  60  m/s. 

5.2  Procedure  Used  to  Measure  Spectra 

Spectra  were  measured  using  the  hot-wire  anemometer  system  and  the  modified  cross¬ 
tunnel  rake  that  were  used  for  the  intensity  measurements.  The  spectral  and  intensity 
experimental  programs  were  combined  so  that  data  for  a  given  spectrum  was  generally 
measured  immediately  before  or  immediately  after  data  for  the  corresponding  intensity. 
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When  measuring  spectra,  each  hot-wire  voltage  was  sampled  65536  times  (i.e.  216  times 
-see  below)  at  a  sampling  frequency  of  20  kHz  and  the  low-pass  filters  were  set  at 
10  kHz.  This  filter  setting  was  used  since  velocity  fluctuations  above  10  kHz  in  the 
free-stream  have  negligible  energy  (Reference  14)  and  the  setting  corresponds  to  the 
maximum  frequency  that  can  be  resolved  when  sampling  at  20  kHz,  due  to  aliasing. 
Spectra  were  obtained  by  processing  hot-wire  voltages,  rather  than  by  processing 
velocities  determined  from  voltages,  as  for  the  intensities,  and  it  was  not  necessary  to 
know  hot-wire  sensitivities.  Power  spectral  densities,  Su,  Sv  and  Sw,  were  computed 
from  the  voltages  using  a  Fast-Fourier-Transform  (FFT)  algorithm  (where  the  number 
of  data  points  required  is  a  power  of  2).  For  each  spectrum,  power  spectral  densities 
were  computed  using  the  65536  (=  216)  consecutively  sampled  voltages,  without  any 
averaging.  For  a  sampling  frequency  of  20  kHz,  the  frequencies  at  which  power 
spectral  densities  were  computed  ranged  from  0.3052  Hz  to  10  kHz  in  increments  of 
0.3052  Hz,  i.e.  20000/65536  Hz,  corresponding  to  32768  (=  2'^)  values  of  frequency. 

Whenever  a  fluctuating  velocity  waveform  is  represented  by  a  finite  number  of  points, 
then  the  mean  velocities  near  the  start  and  finish  of  the  finite  sample  can  be  significantly 
different.  When  computing  spectra  from  these  waveforms,  this  behaviour  can  be 
interpreted  as  a  low-frequency  phenomenon,  coupled  with  high  frequency  effects  due  to 
the  step  change  between  the  beginning  and  end  of  the  sample,  and  it  can  cause  an 
increase  in  spectral  densities  at  low  frequencies.  One  way  of  correcting  for  this 
behaviour  is  to  apply  a  Blackman  window  to  the  data  (Reference  15).  A  Blackman 
window  adjusts  velocities  (or  any  given  measurand)  so  that  their  values  at  the  start  and 
finish  of  the  chosen  range  are  gradually  diminished,  becoming  zero  at  the  extremes. 
The  application  of  a  Blackman  window  also  smooths  the  spectra.  This  is  done  by  first 
smoothing  with  a  Hanning  window  and  then  further  smoothing  with  weights  of  0.16, 
0.68  and  0.16  (Reference  15).  All  spectra  given  in  this  report  have  been  smoothed  with 
a  Blackman  window  and  the  spectra  have  been  normalized  by  dividing  by  the  variance. 


53  Analysis  of  Spectra 

Spectra  corresponding  to  conventional  wind-tunnel-generated  turbulence  should  be 
relatively  smooth,  but  if  the  turbulent  flow  contains  eddies  of  dominant  frequencies, 
such  as  vortices  shed  by  the  fan,  then  these  will  produce  harmonics  or  spikes  on  the 
spectra.  If  there  is  any  significant  “noise”  on  a  fluctuating  turbulent  signal,  caused  by 
say  electrical  interference  or  probe  vibration,  then  this  will  also  cause  spikes.  As  a 
result  of  this  noise,  calculated  turbulence  intensities  will  be  higher  than  actual  intensities 
existing  in  the  tunnel.  By  noting  the  frequency  at  which  the  spikes  occur,  it  is  often 
possible  to  identify  the  source  of  the  irregularities,  e.g.  the  frequency  of  the  spikes  may 
correspond  with  the  frequency  of  the  electrical  mains  supply  or  the  frequency  of 
rotation  of  the  fan.  A  knowledge  of  the  causes  of  the  irregularities  is  valuable  when 
trying  to  eradicate  them. 
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Power  spectral  densities,  Su,  plotted jigainst  frequency,  f,  are  shown  in  Figure  33  at  the 
grid  points  shown  in  Figure  26  for  UN0M  =  30  m/s  for  xT  =  -2.0,  -1.0,  0.0,  and  +1.0  m. 
Corresponding  spectra  for  Sv  and  Sw  are  shown  in  Figures  34  and  35  respectively.  The 
spectra  in  each  of  Figures  33  to  35  are  plotted  together  without  being  identified 
individually.  Figures  36  to  38  show  spectra  for  UNOM  =  60  m/s.  The  spectra  cover  the 
central  region  of  the  flow  in  the  test  section.  The  spectra  show  that  there  are  no  major 
irregularities  in  the  turbulence  signals  over  the  frequency  range  0.3  Hz  to  10  kHz.  For 
all  of  the  u-component  spectra,  there  are  small  but  consistent  spikes  at  about  150  Hz, 
and  the  reason  for  these  is  uncertain.  The  spikes  are  not  associated  with  the  frequency 
of  the  fan  rotation  since  they  are  present  at  about  150  Hz  for  both  UN0M  =  30  and  60  m/s 
and  they  do  not  appear  to  be  caused  by  electrical  interference  from  mains  supply 
(50  Hz).  Similar  spikes  were  found  on  spectra  measured  in  the  LSWT  prior  to  the 
tunnel  being  modified  in  1996,  as  shown  in  Reference  3.  However,  it  was  shown  that 
this  type  of  irregularity  did  not  significantly  affect  measured  values  of  turbulence 
intensities.  There  are  also  some  isolated  spikes  at  the  high-frequency  ends  of  some 
spectra. 


6.  Comparison  of  Flow  Quality  in  the  Original  and  the 
Extended  Test  Sections  of  the  LSWT 


As  indicated  in  Section  1,  the  LSWT  was  modified  in  1996  by  extending  the  length  of 
the  test  section  and  installing  a  new  contraction  (see  Appendix  A  for  details  of  the 
modifications  to  the  tunnel).  This  was  done  to  enable  longer  models  to  be  tested  in  the 
tunnel  and  hopefully  to  obtain  some  improvement  in  the  quality  of  the  flow  in  the  test 
section.  The  flows  in  the  original  and  the  extended  test  sections  are  now  compared  to 
see  how  the  modifications  to  the  tunnel  have  affected  the  flow.  Plots  are  given  for  mean 
velocities,  flow  angles,  turbulence  intensities  and  spectra  for  xT  =  0.0  m  in  the  original 
(see  Reference  3)  and  the  extended  test  sections.  Reference  is  made  to  Table  1,  which 
lists  how  mean  velocities,  flow  angles  and  turbulence  intensities  vary  in  the  central  50% 
of  the  cross-sectional  area  of  the  test  section  between  xT  =  -1.0  and  +1.0  m  for  the 
original  and  the  extended  test  sections.  Further  comparisons  can  be  made  as  required  by 
comparing  other  results  in  this  report  with  corresponding  results  in  Reference  3. 

6.1  Mean  Velocities 

Contours  of  (U  -  U00)/U0  0  for  xT  =  0.0  m  for  UN0M  =  30  and  45  m/s  for  the  extended  and 
the  original  test  section  are  shown  in  Figure  39.  Corresponding  contours  for  UN0M  =  60 
and  75m/s  are  shown  in  Figure  40.  The  distributions  of  corresponding  mean  velocities 
in  the  extended  and  the  original  test  sections  are  very  similar.  Table  1  shows  that  the 
range  of  variation  of  (U  -  U00)/U00  between  xT  =  -1.0  and  +1.0  m  is  generally  slightly 
less  in  the  extended  test  section  over  the  velocity  range.  Overall,  the  modifications  to 
the  tunnel  have  only  very  slightly  improved  the  mean  velocities. 
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Figure  33.  Spectra  for  the  u  component  of  the  turbulence  for  UNOM  =  30  m/s  for 
xt  =  -2-0,  -1 .0,  0.0  and  1 .0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(a)  xT  =  -2.0  and  -1 .0  m. 
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Figure  33  cont’d.  Spectra  for  the  u  component  of  the  turbulence  for  UN0M  =  30  m/s  for 
xT  =  -2.0,  -1 .0,  0.0  and  1 .0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(b)  xT  =  0.0  and  1 .0  m. 
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Figure  34.  Spectra  for  the  v  component  of  the  turbulence  for  UNOM  =  30  m/s  for 
xt  =  _2.0,  -1.0,  0.0  and  1.0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(a)  xT  =  -2.0  and  -1 .0  m. 
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Figure  34  cont’d.  Spectra  for  the  v  component  of  the  turbulence  for  UN0M  =  30  m/s  for 
xT  =  -2.0,  -1.0,  0.0  and  1 .0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(b)  xT  =  0.0  and  1.0  m. 
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Figure  35.  Spectra  for  the  w  component  of  the  turbulence  for  UN0M  =  30  m/s  for 
xt  —  “2.0,  —1 .0,  0.0  and  1 .0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(a)  xr  =  -2.0  and  -1 .0  m. 
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Figure  35  cont’d.  Spectra  for  the  w  component  of  the  turbulence  for  UNOM  =  30  m/s  for 
xT  =  -2.0,  -1.0,  0.0  and  1.0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(b)  xT  =  0.0  and  1 .0  m. 
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Figure  36.  Spectra  for  the  u  component  of  the  turbulence  for  UN0M  =  60  m/s  for 
xt  =  -2-0,  -1 .0,  0.0  and  1 .0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(a)  xT  =  -2.0  and  -1.0  m. 


80 


DSTO-TR-1384 


Power 

Spectral 

Density 

S 


Figure  36  cont’d.  Spectra  for  the  u  component  of  the  turbulence  for  UNOM  =  60  m/s  for 
xT  =  -2.0,  -1.0,  0.0  and  1.0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(b)  xT  =  0.0  and  1.0  m. 
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Figure  37.  Spectra  for  the  v  component  of  the  turbulence  for  UN0M  =  60  m/s  for 
xt  —  ~2.0,  —1.0,  0.0  and  1.0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(a)  xT  =  -2.0  and  -1 .0  m. 
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Figure  37  cont’d.  Spectra  for  the  v  component  of  the  turbulence  for  UNOM  =  60  m/s  for 
xT  =  -2.0,  -1 .0,  0.0  and  1 .0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(b)  xT  =  0.0  and  1 .0  m. 
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Figure  38.  Spectra  for  the  w  component  of  the  turbulence  for  UNOM  =  60  m/s  for 
xT  =  -2.0,  -1 .0,  0.0  and  1 .0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(a)  xT  =  -2.0  and  -1.0  m. 
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Figure  38  cont’d.  Spectra  for  the  w  component  of  the  turbulence  for  UN0M  =  60  m/s  for 
=  -2.0,  -1 .0,  0.0  and  1.0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

tbl  x,  =  0.0  and  1 .0  m. 
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Table  1.  Variation  of  flow  parameters  in  the  original  and  the  extended 
test  sections  of  the  LSWT  for  -1.0  m  <  x.<  +1.0  m. 


Parameter 

^NOM 

LSWT  with  Original 

LSWT  with  Extended 

Test  Section 

Test  Section 

m/s 

Limits  of 

Range  of 

Limits  of 

Range  of 

Variation 

Variation 

Variation 

Variation 

(u-uj/u„ 

30 

-0.6%  to  2.0% 

2.6% 

-1.0%  to  1.4% 

2.4% 

Mean 

Velocity 

Measurements 

(u-uj/u,,. 

45 

-0.8%  to  1 .6% 

2.4% 

-0.6%  to  1 .6% 

2.2% 

(u-uj/u„ 

60 

-0.6%  to  1 .8% 

2.4% 

-0.6%  to  1 .4% 

2.0% 

(u-uj/u„ 

75 

-0.2%  to  2.0% 

2.2% 

-0.8%  to  1 .4% 

2.2% 

30 

-0.8°  to  0.6° 

1.4° 

-0.8°  to  0.4° 

1.2° 

Flow  Angle 

£v 

30 

-0.8°  to  0.4° 

1.2° 

-0.8°  to  0.4° 

1.2° 

Measurements 

60 

-0.4°  to  0.8° 

1.2° 

-0.4°  to  0.8° 

1.2° 

Ev 

60 

-1.0°  to  0.4° 

1.4° 

-1.0°  to  0.4° 

1.4° 

Turbulence 

(5T’/u„ 

30  &  60 

up  to  0.4% 

up  to  0.4% 

Intensity 

Measurements 

30  &  60 

up  to  0.7% 

up  to  0.7% 

30  &  60 

up  to  0.7% 

up  to  0.7% 
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6.2  Flow  Angles 

Contours  of£H  and  £v  for  xT  =  0.0  m  for  UN0M  =  30  m/s  for  the  extended  and  the  original 
test  sections  are  shown  in  Figure  41.  Corresponding  contours  for  UN0M  =  60  m/s  are 
shown  in  Figure  42.  The  distributions  of  corresponding  flow  angles  in  the  extended  and 
the  original  test  sections  are  similar.  Table  1  shows  that  the  range  of  variation  of  flow 
angles  between  xT  =  -1.0  and  +1.0  m  are  approximately  the  same  in  the  original  and  the 
extended  test  sections  for  the  two  velocities  used.  Overall,  the  modifications  to  the 
tunnel  have  had  virtually  no  effect  on  the  flow  angles. 

6.3  Turbulence  Intensities 

Contours  of  (if  )°"7U00,  (v")°7Uon  and  (w^'VU,,,,  for  xT  =  0.0  m  for  UN0M  =  30  m/s  for  the 
extended  and  original  test  sections  are  shown  in  Figure  43.  Corresponding  contours  for 
UN0M  =  60  m/s  are  shown  in  Figure  44.  Table  1  shows  that  the  intensities  between 
xT  =  -1.0  and  +1.0  m  over  the  central  regions  of  the  cross  sections  are  the  same  in  the 
extended  and  the  original  test  sections  at  the  two  velocities.  While  the  distributions  of 
corresponding  intensities  in  the  central  regions  of  the  extended  and  the  original  test 
sections  are  similar,  the  gradients  of  the  intensities  and  absolute  values  of  intensities 
near  the  walls  of  the  test  sections  are  often  significantly  lower  for  the  extended  test 
section.  This  could  be  due  to  the  fact  that  the  flow  near  the  contraction  walls  has  had  to 
negotiate  different  curved  paths  in  the  two  cases  because  the  shapes  of  the  two 
contractions  are  different.  However,  the  reasons  for  the  observed  flow  behaviour  in  the 
current  investigation  are  probably  quite  complex,  since  the  two  contractions  have 
different  concave  and  also  different  convex  radii  of  curvature,  and  it  is  not  possible  to 
give  a  precise  explanation  for  the  changes. 

6.4  Spectra 

Spectra  for  the  u-,  v-  and  w-components  of  the  turbulence  for  xT  =  0.0  m  for 
tf  „  =  30  m/s  for  the  extended  and  the  original  test  sections  are  shown  in  Figure  45. 
Corresponding  spectra  for  UN0M  =  60  m/s  are  shown  in  Figure  46.  Corresponding 
spectra  in  the  extended  and  the  original  test  sections  are  virtually  the  same,  indicating 
that  the  modifications  to  the  tunnel  have  had  little  effect  on  the  spectra. 

6.5  Overall  Effect  of  Modifying  the  Contraction  and  the  Length  of  the  Test 
Section  of  the  LSWT  in  1996 

The  installation  of  the  extended  test  section  and  the  new  contraction  has  had  only  a 
minimal  effect  on  improving  the  quality  of  the  flow  in  the  test  section.  However,  since 
the  contraction  ratio  for  the  new  contraction  was  the  same  as  that  for  the  original 
contraction,  i.e.  4:1,  it  was  not  expected  that  simply  altering  the  shape  of  the  contraction 
would  result  in  a  large  improvement  in  the  quality  of  the  flow  in  the  test  section.  The 
main  reason  for  modifying  the  tunnel  in  1996  was  to  extend  the  length  of  the  test 
section  to  enable  longer  models  to  be  tested. 
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Figure  40  cont’d.  Contours  of  velocity  deviation  for  UNOM  =  60  and  75  m/s  for  xT 
(b)  Data  taken  in  the  LSWT  prior  to  modifications  to  the  tunnel  in  1996. 
The  flow  is  out  of  the  page. 


Figure  43  cont'd.  Contours  of  u-  v-  and  w-component  intensities  for  UNOM  =  30  m/s  for  x 
(b)  Data  taken  in  the  LSWT  prior  to  modifications  to  the  tunnel  in  1996. 
The  flow  is  out  of  the  page. 
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Figure  45.  Spectra  for  the  u,  v  and  w  components  of  the  turbulence  for  UN0M  =  30  m/s 
for  xT  =  0.0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(a)  Data  taken  in  the  LSWT  with  the  extended  test  section  and  new  contraction. 


Figure  45  cont’d.  Spectra  for  the  u,  v  and  w  components  of  the  turbulence  for  UNOM  =  30  m/s 
for  xT  =  0.0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(b)  Data  taken  in  the  LSWT  prior  to  modifications  to  the  tunnel  in  1996. 
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Figure  46.  Spectra  for  the  u,  v  and  w  components  of  the  turbulence  for  UN()M  =  60  m/s 
for  xT  =  0.0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(a)  Data  taken  in  the  LSWT  with  the  extended  test  section  and  new  contraction. 


Figure  46  cont’d.  Spectra  for  the  u,  v  and  w  components  of  the  turbulence  for  UNOM  =  60  m/s 
for  xT  =  0.0  m  for  yT  and  zT  locations  denoted  by  filled  circles  in  Figure  26. 

(b)  Data  taken  in  the  LSWT  prior  to  modifications  to  the  tunnel  in  1996. 


103 


DSTO-TR-1384 


7.  Assessment  of  Flow  Quality  in  the  LS  WT 

It  has  been  shown  that,  like  nearly  all  fully-enclosed  flows,  the  flow  in  the  extended  test 
section  of  the  LSWT  has  some  non-uniformity  in  terms  of  mean-velocity  deviations, 
flow  angularity  and  turbulence  intensities.  In  the  following,  flow  non-uniformities  for 
the  tunnel  are  quoted  for  approximately  the  central  50%  of  the  cross-sectional  area  of 
the  test  section. 

7.1  Mean  Velocities 

F°r  UNnM  =  30  m/s,  between  xT  =  -2.0  and  +1 .0  m,  longitudinal  mean  velocities  deviated 
between  -1.0%  and  +1.4%  from  their  values  at  the  centre  of  the  test  section  (xT  =  0.0, 
yT  =  0.0  and  z^  =  0.0  m)  (Figure  8).  Corresponding  deviations  for  other  values  of  UN0M 
were  from  -1.2%  to  +1.6%  for  UNOM  =  45  m/s  (Figure_9),  from  -1.0%  to  +1.4%  for 
Unom  =  60  m/s  (Figure  1 0)  and  from  - 1 .2%  to  +1.4%  for  UN0M  =  75  m/s  (Figure  1 1 ). 

In  Reference  10  it  is  indicated  that  acceptable  longitudinal  mean-velocity  deviations 
across  the  test  section  are  often  quoted  in  the  range  ±0.20%  to  ±0.30%  variation  from 
the  average  value  of  the  mean  velocity  (not  stated  where  in  the  test  section). 
Reference  16  indicates  that  acceptable  mean-velocity  variations  across  the  test  section 
of  a  high  performance  tunnel  are  less  than  ±0.2%.  In  the  8  m  x  6  m  test  section  of  the 
German-Dutch  large  low-speed  wind  tunnel,  the  total-pressure  deviations  (expressed  as 
pressure  coefficients)  are  less  than  0.3%  at  90  m/s,  and  the  static-pressure  deviations 
(expressed  as  pressure  coefficients)  are  well  within  0.6%  at  100  m/s  (Reference  14). 
The  deviations  in  the  total  and  static  pressures  are  relative  to  the  total  and  static 
pressures  respectively  at  the  centre  of  the  tunnel  test  section.  (Unfortunately,  it  has  not 
been  possible  to  determine  variations  in  mean  velocity  from  the  published  data  since  the 
given  deviations  in  total  pressure  apply  to  a  region  and  not  specific  locations,  so  that 
total  pressures  and  static  pressures  cannot  be  matched).  According  to  Reference  17, 
velocities  should  not  deviate  from  the  mean  value  by  more  than  ±0.75%.  The  five 
metre  low-speed  wind  tunnel  at  the  Defence  Evaluation  and  Research  Agency  (formerly 
Royal  Aircraft  Establishment)  has  mean-velocity  variations  of  about  ±0.12%  although 
the  design  aim  was  only  ±0.05%  (Reference  18). 

12  Flow  Angles 

For  UNOM  =  30  m/s,  between  xT  =  -2.0  and  +1.0  m,  horizontal  flow  angles  varied  from 
-0.8°  to  +0.4°  (Figure  21)  and  vertical  flow  angles  varied  from  -1.0°  to  +0.4° 
(Figure  22).  For  UNOM  =  60  m/s,  horizontal  flow  angles  varied  from  -0.6°  to  +0.8° 
(Figure  23)  and  vertical  flow  angles  varied  from  -1.2°  to  +0.4°  (Figure  24). 

According  to  Reference  10,  it  would  be  desirable  to  have  horizontal  and  vertical  flow 
angle  variations  of  less  than  ±0.10°,  but  if  this  cannot  be  achieved,  then  the  variations 
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should  be  held  to  within  ±0.20°.  In  Reference  16  it  is  indicated  that  acceptable  flow 
angle  variations  for  a  high  performance  tunnel  are  ±0.1°.  For  the  German-Dutch  wind 
tunnel  (Reference  14),  the  horizontal  and  vertical  flow  angles  vary  by  less  than  0.2° 
across  the  central  region  of  the  8  m  x  6  m  test  section  for  a  reference  free-stream 
velocity  of  90  m/s.  According  to  Reference  17,  horizontal  and  vertical  flow  angles 
should  not  deviate  from  the  axial  direction  by  more  than  ±0.25°. 

13  Turbulence  Intensities  and  Spectra 

For  UNnM  =  30  and  60  m/s,  for  xT  =  -2.0,  -1.0,  0.0  and  +1 .0  m,  values  of  (u2)°'VU  0  were 
generally  below  about  0.4%,  whereas  corresponding  values  of  (v2)  /U00  and  (w2)  /U00 
were  generally  below  about  0.7%  (Figures  27  to  32).  Spectra  for  the  u-  v-  and  w- 
components  of  the  turbulence  were  relatively  smooth  over  the  frequency  range  0.3  Hz 
to  10  kHz,  with  no  significant  spikes  associated  with  irregularities  in  the  flow 
(Figures  33  to  38). 

Reference  13  indicates  that  a  root-mean-square  u-component  fluctuation  of  0.1%  of  the 
mean  velocity  is  often  quoted  as  low  enough  for  most  experiments  and  that  the  best 
low-turbulence  research  tunnels  have  a  turbulence  intensity  of  the  order  of  0.02%  at  low 
speeds.  In  Reference  10  it  is  indicated  that  tunnels  used  for  developmental  testing  can 
have  longitudinal  turbulence  intensities  as  high  as  0.5%,  and  that  tunnels  used  for 
research  in  boundary  layers  and  boundary  layer  transition  should  have  lateral  values  of 
intensity  below  about  0.05%.  In  the  German-Dutch  wind  tunnel  (Reference  14),  the 
longitudinal  intensities  are  quoted  to  be  less  than  about  0.08%  and  the  lateral  intensities 
less  than  about  0.15%  in  the  central  region  of  the  8  m  x  6  m  test  section  for  a  reference 
free-stream  velocity  of  75  m/s.  However,  for  a  velocity  of  about  40  m/s,  the 
longitudinal  intensities  are  0.02%,  which  is  comparable  with  the  best  low  turbulence 
wind  tunnels. 

7.4  Adequacy  of  Flow  Quality  in  the  LSWT 

The  foregoing  suggests  that  although  the  quality  of  the  flow  in  the  extended  test  section 
of  the  LSWT  is  acceptable,  it  could  be  improved.  However,  there  is  no  point  spending  a 
large  amount  of  money  on  the  tunnel  to  improve  the  flow  if  it  is  adequate  for  the 
investigations  likely  to  be  undertaken.  The  question  of  whether  the  quality  of  the  flow 
is  adequate  depends  upon  the  types  of  tests  likely  to  be  undertaken  in  the  tunnel. 

For  flow-visualization  investigations,  where  qualitative  rather  than  quantitative  answers 
are  often  required,  the  quality  of  the  flow  is  not  all  that  important  and  the  tunnel  should 
be  adequate  for  this  type  of  work. 

For  experiments  in  which  boundary-layer  transition  phenomena  are  important,  the  free- 
stream  turbulence  intensities  should  be  low  since  high  intensities  can  cause  early 
transition  and  this  can  lead  to  difficulties  in  interpreting  experimental  results. 
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Longitudinal  turbulence  intensities  of  0.1%  to  0.5%,  similar  to  those  in  the  LSWT,  are 
classed  as  moderate  turbulence  levels  (Reference  19),  so  that  some  lack  of  similarity 
may  result  if  the  tunnel  was  used  for  these  types  of  experiments. 

Accurate  force,  moment  and  surface  pressure  measurements  taken  on  a  model  of  say  an 
aircraft  are  needed  in  a  number  of  areas,  including  to  determine  loading  programs  for 
fatigue  tests  on  full-scale  aircraft,  to  provide  data  for  the  development  of  flight-dynamic 
models  of  aircraft  behaviour,  and  to  provide  data  for  the  verification  of  computational- 
fluid-dynamic  (CFD)  codes.  If  there  are  imperfections  in  the  flow,  then  this  may  lead  to 
data  being  produced  that  does  not  represent  the  full-scale  aircraft  to  the  accuracy 
required.  For  example,  if  the  angle  of  the  flow  incident  on  one  of  the  wings  of  a  model 
of  a  typical  fighter  aircraft  is  1  °  different  from  that  on  the  other  wing,  due  to  the  flow 
asymmetry  in  the  test  section,  then  potentially  this  could  produce  an  error  of  the  order 
of  5  to  10%  in  the  maximum  measured  rolling  moment  coefficient  for  the  aircraft. 
However,  to  provide  accurate  data  from  wind-tunnel  tests  to  enable  the  performance  of 
a  full-scale  aircraft  to  be  predicted,  irregularities  in  the  flow  is  only  one  of  the  factors 
that  needs  to  be  considered.  Reynolds  number  effects  and  the  size  of  the  model  with 
respect  to  the  size  of  the  test  section  are  also  important. 

Overall,  the  quality  of  the  flow  in  the  LSWT  is  considered  to  be  adequate  to  enable  tests 
to  be  carried  out  with  sufficient  accuracy  to  meet  the  requirements  of  the  majority  of  the 
test  programs  undertaken  in  the  tunnel. 


8.  Possible  Causes  of  Flow  Non  Uniformities  and  Possible 
Methods  for  Improving  the  Flow 

Flow  irregularities  in  the  test  sections  of  wind  tunnels  may  arise  from  a  number  of 
different  flow  mechanisms  and  generally  it  is  not  possible  to  identify  with  certainty  the 
causes  of  these  irregularities  simply  by  analysing  measured  calibration  data.  Variations 
in  mean  velocity  can  be  caused  by  flow  separations  in  the  tunnel  return  circuit,  flow 
separations  in  the  contraction,  poor  corner  vane  design  or  incorrectly  set  vanes 
(resulting  in  the  vanes  over  or  under  turning  the  flow),  and  poor  design  of  the  fan  or 
straightening  vanes  (resulting  in  a  rotation  of  the  whole  flow  downstream  of  the  fan). 
Turbulence  can  be  generated  in  the  boundary  layers  on  the  tunnel  walls  and  from  the 
wakes  of  the  fan  blades,  the  turning  vanes  and  the  honeycombs.  Angularity  in  the  flow 
can  persist  into  the  test  section  if  the  honeycomb  has  cells  of  an  inappropriate  size  and 
aspect  ratio  (see  below).  Despite  the  uncertainties  involved  in  making  a  precise  assess¬ 
ment  of  the  causes  of  flow  irregularities  in  a  tunnel,  it  is  possible  to  identify  some  of  the 
likely  causes  for  the  irregularities  in  the  extended  test  section  of  the  LSWT  at  PSL. 

If  the  contraction  in  a  wind  tunnel  is  not  well  designed  then  adverse  pressure  gradients 
can  occur  near  the  inlet  and  this  can  lead  to  either  continuous  flow  separation,  or 
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intermittent  flow  separation,  depending  on  the  magnitude  of  the  adverse  pressure 
gradient  (Reference  10).  This  separation  can  have  an  adverse  effect  on  the  quality  of 
the  flow  in  the  test  section.  Recent  wool-tuft  studies  near  the  inlet  of  the  new 
contraction  did  not  show  any  evidence  of  flow  separation.  This  suggests  that  the  small 
flow  irregularities  that  occur  in  the  extended  test  section  must  come  from  other  sources 
in  the  tunnel  flow  circuit. 

Honeycombs  and  screens  are  known  to  improve  the  quality  of  the  flow  in  the  test 
section  of  a  tunnel,  both  in  terms  of  improved  longitudinal  and  lateral  mean-velocity 
distributions  and  reduced  longitudinal  and  lateral  turbulence  intensities.  Screens  reduce 
longitudinal  components  of  mean  velocities  and  intensities  more  than  lateral 
components,  as  for  a  contraction,  whereas  honeycombs  reduce  lateral  components  of 
mean  velocities  and  intensities  more  than  longitudinal  components.  The  honeycomb  in 
the  LSWT  has  triangular  cells  whose  frontal  dimensions  are  48,  41  and  41  mm  and  the 
depth  of  the  cells  is  127  mm  (approximate  dimensions),  as  shown  in  Figure  47.  The 
diameter  of  a  circle  having  the  same  cross-sectional  area  of  a  cell  is  about  32  mm  so  that 
the  ratio  of  the  depth  of  the  cells  to  the  effective  diameter  of  the  cells  is  about  4:1, 
which  is  somewhat  lower  than  the  recommended  aspect  ratio  of  about  6:1  to  8:1 
(References  10  and  13).  The  sizes  of  the  cells  in  the  honeycomb  are  also  larger  than 
normal.  In  Reference  20,  a  study  was  made  of  honeycombs  in  wind  tunnels  and  it  was 
found  that  the  honeycombs  had  cell  sizes  ranging  from  1.6  mm  (1/16  inch)  to  9.5  mm 
(3/8  inch)  and  aspect  ratios  varying  between  6:1  and  8:1.  The  existing  honeycomb  is 
not  continuous  as  it  has  been  installed  in  the  tunnel  in  16  sections.  Between  each 
section  there  is  a  solid  edge  about  15  mm  thick  and  these  edges  may  have  a  detrimental 
effect  on  the  flow  in  the  test  section  (Reference  10).  If  it  is  decided  that  a  new 
honeycomb  is  needed  (see  below),  then  it  should  have  cells  with  a  smaller  frontal  area 
and  a  larger  aspect  ratio,  and  follow  the  design  dimension  guidelines  given  above.  It 
would  be  preferable  to  use  a  honeycomb  having  hexagonal  cells  since  the  pressure  loss 
for  this  type  of  honeycomb  is  less  than  for  honeycombs  having  square  or  circular  cross 
sections  (Reference  10).  It  would  probably  be  necessary  to  use  some  type  of  supporting 
structure  to  prevent  the  honeycomb  deflecting,  and  it  should  be  less  intrusive  than  the 
rigid  supporting  structure  presently  in  the  tunnel.  Currently  the  tunnel  does  not  have 
any  screens.  If  it  is  decided  to  install  them,  an  open-to-total  area  ratio  of  the  order  of 
57%  is  preferable  (Reference  21). 

Most  likely  there  is  not  a  single  cause  of  the  flow  non  uniformities  and  it  may  be 
necessary  to  change  more  than  one  component  of  the  tunnel  to  obtain  a  large  improve¬ 
ment  in  the  quality  of  the  flow  in  the  test  section.  It  is  preferable  that  the  modifications 
to  the  tunnel  be  carried  out  in  stages.  Extending  the  test  section  and  installing  a 
redesigned  contraction  has  had  only  a  minimal  effect  on  improving  the  flow.  It  will  be 
necessary  to  make  additional  modifications  to  the  tunnel  to  improve  the  flow  further.  If 
it  is  decided  to  install  a  new  honeycomb,  then  preformed  aluminium  hexagonal 
honeycomb  or  stacks  of  thin-walled  hexagonal  metal  tubing  could  be  inserted  into  the 
existing  honeycomb  structure  to  replace  the  current  triangular  cells  and  minimize  cost. 
Modifying  the  honeycomb  and/or  installing  screens  would  increase  losses  and  reduce 
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(a)  view  of  complete  honeycomb,  showing  details  of  the  16  sections 


(b)  close-up  view  of  the  honeycomb,  showing  frontal  details  of  the  triangular  cells 
Figure  47.  Honeycomb  used  in  the  low-speed  wind  tunnel. 
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the  maximum  velocity  obtainable  in  the  test  section.  Additional  cleaning  would  also  be 
necessary  as  both  the  modified  honeycomb  and  screens  could  accumulate  more  dust 
which  could  change  the  quality  of  the  flow  with  time  (Reference  10).  A  new 
honeycomb  may  provide  sufficient  improvement  in  the  flow  without  screens.  Further 
studies  would  be  needed  to  determine  the  best  way  to  proceed. 


9.  Concluding  Remarks 


A  detailed  calibration  of  the  flow  in  the  extended  test  section  of  the  low-speed  wind 
tunnel  is  given  in  this  report.  An  extensive  range  of  longitudinal  mean-flow  velocities, 
flow  angles,  turbulence  intensities  and  turbulence  spectra  were  measured  at  up  to  155 
grid  points  across  the  flow  for  longitudinal  positions,  xT,  of  -2.0,  -1.0,  0.0  and  +1.0  m 
from  the  centre  of  the  turntables  in  the  test  section  and  for  nominal  free-stream 
velocities,  UNOM,  ranging  from  30  to  75  m/s. 

Longitudinal  mean-flow  velocities  were  measured  using  pressure  probes  fitted  to  a  rake 
which  spanned  the  test  section.  Flow  angles  were  measured  using  a  six-hole  probe 
mounted  on  a  specially-designed  traversing  mechanism.  Turbulence  intensities  and 
power  spectral  densities  were  obtained  using  constant-temperature  hot-wire 
anemometers. 

Considering  approximately  the  central  50%  of  the  cross-sectional  area  of  the  test 
section,  between  a  length  of  xT  =  -2.0  and  +1.0  m,  the  following  flow  non-uniformities 
were  observed. 

•  Mean-Flow  Velocities 

For  UNOM  =  30  m/s,  longitudinal  mean  velocities  deviated  between  -1.0%  and  +1.4% 
from  the  mean  value  at  the  centre  of  the  test  section  (xT  =  0.0,  yT  =  0.0  and  zT  =  0.0  m). 
Corresponding  deviations  for  other  values  of  UN0M  were  from  -1.2%  to  +1.6%  for 
UN0M  =  45  m/s,  from  -1.0%  to  +1 .4%  for  UN0M  =  60  m/s,  and  from  -1.2%  to  +1.4%  for 
UN0M  =  75  m/s.  No  significant  periodic  variations  in  mean  velocities  were  found. 

•  Flow  angles 

For  UNOM  =  30  m/s,  horizontal  flow  anglesjvaried  from  -0.8°  to  +0.4°  and  vertical  flow 
angles  varied  from  -1.0°  to  +0.4°.  For  UN0M  =  60  m/s,  horizontal  flow  angles  varied 
from  -0.6°  to  +0.8°  and  vertical  flow  angles  varied  from  -1.2°  to  +0.4°. 

•  Turbulence  Intensities  and  Spectra 

For  UNOM  =  30  and  60  m/s,  at  xT  =  -2.0,  -1.0,  0.0  and  +1.0  m,  the  root-mean-square 
u-component  intensities  were  generally  below  about  0.4%,  whereas  corresponding  v- 
and  w-component  intensities  were  generally  below  about  0.7%.  Spectra  for  the  u-  v- 
and  w-components  of  the  turbulence  were  relatively  smooth  over  the  frequency  range 
0.3  Hz  to  10  kHz,  with  no  significant  spikes  associated  with  irregularities  in  the  flow. 
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The  installation  of  the  extended  test  section  and  new  contraction  in  1996  has  slightly 
improved  the  quality  of  the  flow  in  the  test  section,  but  some  variations  in  mean 
velocities,  flow  angles,  turbulence  intensities  and  spectra  still  remain.  From  the  present 
calibration  measurements  it  is  not  possible  to  identify  with  certainty  the  causes  of  these 
irregularities.  Wool-tuft  studies  near  the  inlet  to  the  new  contraction  did  not  show  any 
evidence  of  flow  separation,  and  it  is  unlikely  that  the  flow  irregularities  result  from  any 
flow  separation  in  this  area.  The  honeycomb  has  larger  cells  with  a  lower  aspect  ratio 
(cell  depth  divided  by  effective  cell  diameter)  than  recommended  for  this  type  of  tunnel, 
and  this  could  contribute  to  the  irregularities.  If  further  improvement  of  flow  quality  is 
required,  the  honeycomb  could  be  replaced  with  one  having  smaller  cells  that  have  a 
larger  aspect  ratio.  The  tunnel  does  not  have  any  screens  and  it  may  be  necessary  to 
install  them  to  reduce  turbulence  intensities.  However,  this  will  result  in  a  reduction  in 
the  maximum  velocity  obtainable  in  the  test  section  and  increased  maintenance, 
particularly  cleaning. 
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Appendix  A:  Modifications  to  the  LSWT  in  1996 


To  enable  longer  models  to  be  tested,  an  extended  test  section  and  a  new  contraction 
were  installed  in  the  LSWT  in  1996.  It  was  anticipated  that  these  modifications  would 
also  produce  a  slight  improvement  in  the  quality  of  the  flow  in  the  test  section.  The 
tunnel  was  modified  according  to  the  design  recommendations  given  in  1986  by 
Watmuff  (see  Reference  2). 

A.l  Details  of  the  Extended  Test  Section 

The  length  of  the  test  section  was  increased  from  4.12  m  to  6.56  m  by  installing  a  liner 
in  the  tunnel  upstream  of  the  original  test  section.  The  cross-sectional  dimensions  of  the 
extended  test  section  were  the  same  as  those  for  the  existing  test  section. 


A.2  Details  of  the  New  Contraction 

The  cross-sectional  dimensions  of  the  new  contraction  at  its  inlet  and  outlet  were  the 
same  as  those  for  the  original  contraction. 

The  procedure  used  by  Watmuff  to  select  a  new  contraction  shape  for  the  LSWT  is 
described  in  detail  in  Reference  2,  and  summarized  in  the  following. 

Watmuff  indicated  that  a  variety  of  axi-symmetric  contraction  shapes  can  be  specified 
by  two  simple  polynomial  arcs  which  are  joined  together  at  an  inflection  point  (see 
Reference  2),  as  shown  in  Figure  Al.  Using  the  notation  shown  in  this  figure,  the 
equations  for  the  two  arcs,  AB  and  BC,  are  respectively 

R=Ri-A1x"  (Al) 

and 

R=  Rq  +  A2(Lc-x)n  (A2) 


where  R  is  the  distance  from  the  tunnel  centreline  to  the  surface  of  an  axi-symmetric 
contraction  and  x  is  the  distance  from  the  contraction  inlet  in  the  longitudinal  direction. 
Subscripts  i  and  o  refer  to  the  inlet  and  the  outlet  of  the  contraction  respectively.  Lc 
denotes  the  contraction  length  and  n  is  an  exponent  which  controls  the  maximum 
steepness  of  the  contraction  shape.  By  matching  the  gradient  of  the  arcs  at  the 
inflection  point,  the  constants  A]  and  A2  are  given  by 
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Figure  A1 .  Notation  for  contraction  wall  contour  constructed  from  two  arcs 
(based  on  a  figure  given  in  Reference  2). 


and 


4  R:  -  R 

A,  =  — ! - o- 


(A3) 


A2  - 


R.  -  R„ 


(Lc-Xjyn-')Lc 


(A4) 


respectively.  X{  is  the  longitudinal  distance  from  the  contraction  inlet  to  the  inflection 
point.  Different  contraction  shapes  are  obtained  using  different  values  of  L  n  R  R 
and  X;.  1,0 

A  question  that  had  to  be  answered  by  Watmuff  was  what  values  of  these  parameters 
were  the  “best”  to  use  for  the  new  contraction  in  the  LSWT.  He  decided  that  the  shape 
of  the  new  contraction  should  be  selected  on  the  basis  of  the  likely  pressure  distribution 
along  its  length.  He  indicated  that  a  reasonable  criterion  for  the  selection  of  a  new 
contraction  shape  is  to  minimize  any  adverse  pressure  gradient  near  the  contraction  inlet 
while  not  significantly  increasing  the  adverse  pressure  gradient  near  the  contraction 
outlet,  compared  with  the  then  existing  contraction.  It  is  known  that  adverse  pressure 
gradients  at  the  inlet  to  a  contraction  should  be  as  small  as  possible  to  reduce  the 
likelihood  of  flow  separation  at  the  contraction  inlet,  which  could  adversely  affect  the 
flow  quality  in  the  test  section  (Reference  10).  Watmuff  developed  a  numerical  method 
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for  the  solution  of  Stoke’ s  stream  function,  enabling  pressure  distributions  for  idealized 
axi-symmetric  flow  to  be  predicted.  He  used  the  technique  to  predict  idealized  flow 
through  the  then  existing  LSWT  contraction  shape,  using  the  approximation  of  axial 
symmetry  for  the  octagonal  cross  section  of  the  contraction,  and  also  neglecting  viscous 
forces.  Calculated  theoretical  pressure  distributions  were  compared  with  experimental 
measurements  obtained  in  a  pilot  tunnel  using  a  scale  model  of  the  contraction  specially 
built  to  measure  the  pressure  distributions.  Reasonable  agreement  was  obtained,  which 
meant  that  the  method  could  be  used  to  determine  a  new  contraction  shape  for  the 
LSWT. 

Watmuff  calculated  the  likely  pressure  distributions  along  a  range  of  axi-symmetric 
contractions  having  different  values  of  n,  Xj/Lc  and  Lc/D;,  where  D;  is  the  inlet  diameter 
of  an  axi-symmetric  contraction.  He  showed  that  for  n  =  3.0,  Xj/Lc  =  0.5  and 
Lc/D;  =  1.0,  the  calculated  pressure  gradients  at  the  inlet  and  outlet  of  a  contraction 
were  about  20%  and  90%  respectively  of  those  for  the  then  existing  contraction  Based 
on  this  finding,  the  dimensions  of  the  eight  surfaces  of  a  new  octagonal  contraction  for 
the  LSWT  were  made  to  correspond  to  n  =  3.0,  Xj/Lc  =  0.5  and  Lc/Dj  =  1.0,  where 
D  j  =  2R;  =  5187.73  mm  is  the  diameter  of  a  circle  having  the  same  cross-sectional  area 
as  the  inlet  to  the  contraction  in  the  LSWT.  Similarly,  D0  =  2R0  =  2593.9  mm  is  the 
diameter  of  a  circle  having  the  same  cross-sectional  area  as  the  outlet  to  the  contraction 
in  the  LSWT 

A  diagrammatic  view  of  the  new  and  original  (prior  to  being  modified  in  1996) 
contractions  is  shown  in  Figure  A2.  The  curved  surfaces  of  the  new  contraction, 
labeled  as  curves  1,  2,  3  and  4  in  Figure  A2,  are  defined  in  terms  of  mathematical 
relationships.  Curves  2  and  4  for  the  comer  fillets  are  the  same.  It  is  convenient  to  use 
a  different  coordinate  system  for  each  type  of  curved  surface  when  specifying 
dimensions.  The  coordinate  system  for  curve  1  is  shown  by  thick  lines  in  Figure  A2. 
The  coordinate  systems  for  curves  2,  3  and  4  have  the  same  origin  and  the  same  x  axis 
as  the  coordinate  system  for  curve  1,  but  are  rolled  about  the  x  axis  by  45°,  90°  and 
135°  respectively,  as  shown  in  Figure  A2  (z  axes  not  shown  to  avoid  confusion). 

Parameters  have  the  following  values  for  the  different  curves: 

•  Curve  1  (upper  and  lower  surfaces) 

Lc  =  5187.7  mm 
n  =  3.0 

Rj  =  2133.6  mm 
R0  =  1066.8  mm 
Xj  =  2593.9  mm 


3  In  this  Appendix,  dimensions  are  given  to  an  accuracy  of  0. 1  mm  since  it  may  be  necessary  in  the  future 
to  set  up  the  idealized  boundaries  of  the  flow  when  predicting  the  flow  through  the  tunnel  using 
computational-fluid-dynamic  (CFD)  techniques. 
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Figure  A2.  Diagrammatic  representation  of  the  new  and  original  contractions,  showing  the 
coordinate  systems  used  to  define  the  curved  surfaces  on  the  new  contraction. 
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•  Curves  2  and  4  (surfaces  forming  the  comer  fillets) 

Lc  =  5187.7  mm 

n  =  3.0 

Rj  =  2694.1  mm 
R0  =  1347.0  mm 
Xj  =2593.9  mm 

•  Curve  3  (side  surfaces) 

Lr  =  5187.7  mm 

n  =  3.0 

Rj  =  2743.2  mm 
R0  =  1371.6  mm 
Xj  =  2593.9  mm 

Using  these  parameters,  x  and  R  coordinates  corresponding  to  the  four  different  curves 
were  computed  using  equations  A1  to  A4  and  the  coordinates  are  given  in  Tables  A1  to 
A4.  Coordinates  were  computed  for  x  increments  of  50  mm.  Positive  and  negative 
values  of  z  defining  points  on  a  surface  where  it  intersects  with  adjoining  surfaces  are 
also  shown  in  Tables  A1  to  A4.  All  dimensions  given  in  Tables  A1  to  A4  apply  to  the 
inner-flow  surfaces  of  the  contraction.  The  contraction  was  constructed  to  these 
coordinates. 
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Table  A1.  Dimensions  of  the  Upper  and  Lower  Surfaces  of  the  Contraction  (Curve  1 ). 


X 

mm 

R 

mm 

z(+) 

mm 

z(-) 

mm 

X 

mm 

R 

mm 

z(+) 

mm 

Z(-) 

mm 

0.0 

2133.6 

1676.4 

-1676.4 

2600.0 

1596.4 

1254.3 

-1254.3 

50.0 

2133.6 

1676.4 

-1676.4 

2650.0 

1566.3 

1230.7 

-1230.7 

100.0 

2133.6 

1676.4 

-1676.4 

2700.0 

1537.4 

1207.9 

-1207.9 

150.0 

2133.5 

1676.3 

-1676.3 

2750.0 

1509.6 

1186.1 

-1186.1 

200.0 

2133.4 

1676.2 

-1676.2 

2800.0 

1482.9 

1165.1 

-1165.1 

250.0 

2133.1 

1676.0 

-1676.0 

2850.0 

1457.3 

1145.0 

-1145.0 

300.0 

2132.8 

1675.8 

-1675.8 

2900.0 

1432.8 

1125.7 

-1125.7 

350.0 

2132.3 

1675.4 

-1675.4 

2950.0 

1409.3 

1107.3 

-1107.3 

400.0 

2131.6 

1674.9 

-1674.9 

3000.0 

1386.8 

1089.7 

-1089.7 

450.0 

2130.8 

1674.2 

-1674.2 

3050.0 

1365.4 

1072.8 

-1072.8 

500.0 

2129.8 

1673.4 

-167.3.4 

3100.0 

1344.9 

1056.7 

-1056.7 

550.0 

2128.5 

1672.4 

-1672.4 

3150.0 

1325.4 

1041.4 

-1041.4 

600.0 

2127.0 

1671.2 

-1671.2 

3200.0 

1306.8 

1026.8 

-1026.8 

650.0 

2125.2 

1669.8 

-1669.8 

3250.0 

1289.2 

1012.9 

-1012.9 

700.0 

2123.1 

1668.2 

-1668.2 

3300.0 

1272.4 

999.7 

-999.7 

750.0 

2120.7 

1666.3 

-1666.3 

3350.0 

1256.5 

987.2 

-987.2 

800.0 

2118.0 

1664.1 

-1664.1 

3400.0 

1241.4 

975.4 

-975.4 

850.0 

2114.8 

1661.7 

-1661.7 

3450.0 

1227.2 

964.2 

-964.2 

900.0 

2111.3 

1658.9 

-1658.9 

3500.0 

1213.7 

953.6 

-953.6 

950.0 

2107.4 

1655.8 

-1655.8 

3550.0 

1201.1 

943.7 

-943.7 

1000.0 

2103.0 

1652.4 

-1652.4 

3600.0 

1189.1 

934.3 

-934.3 

1050.0 

2098.2 

1648.6 

-1648.6 

3650.0 

1177.9 

925.5 

-925.5 

1100.0 

2092.9 

1644.4 

-1644.4 

3700.0 

1 167.4 

917.3 

-917.3 

1150.0 

2087.1 

1639.9 

-1639.9 

3750.0 

1157.6 

909.6 

-909.6 

1200.0 

2080.8 

1634.9 

-1634.9 

3800.0 

1148.5 

902.4 

-902.4 

1250.0 

2073.9 

1629.5 

-1629.5 

3850.0 

1140.0 

895.7 

-895.7 

1300.0 

2066.5 

1623.6 

-1623.6 

3900.0 

1132.1 

889.5 

-889.5 

1350.0 

2058.4 

1617.3 

-1617.3 

3950.0 

1124.8 

883.7 

-883.7 

1400.0 

2049.7 

1610.5 

-1610.5 

4000.0 

1118.0 

878.4 

-878.4 

1450.0 

2040.4 

1603.2 

-1603.2 

4050.0 

1111.8 

873.6 

-873.6 

1500.0 

2030.4 

1595.4 

-1595.4 

4100.0 

1106.1 

869.1 

-869.1 

1550.0 

2019.8 

1587.0 

-1587.0 

4150.0 

1101.0 

865.0 

-865.0 

1600.0 

2008.4 

1578.0 

-1578.0 

4200.0 

1096.3 

861.3 

-861.3 

1650.0 

1996.3 

1568.5 

-1568.5 

4250.0 

1092.0 

858.0 

-858.0 

1700.0 

1983.4 

1558.4 

-1558.4 

4300.0 

1088.2 

855.0 

-855.0 

1750.0 

1969.8 

1547.7 

-1547.7 

4350.0 

1084.8 

852.3 

-852.3 

1800.0 

1955.4 

1536.3 

-1536.3 

4400.0 

1081.7 

849.9 

-849.9 

1850.0 

1940.1 

1524.4 

-1524.4 

4450.0 

1079.1 

847.8 

-847.8 

1900.0 

1924.0 

1511.7 

-1511.7 

4500.0 

1076.7 

846.0 

-846.0 

1950.0 

1907.0 

1498.3 

-1498.3 

4550.0 

1074.7 

844.4 

-844.4 

2000.0 

1889.1 

1484.3 

-1484.3 

4600.0 

1073.0 

843.1 

-843.1 

2050.0 

1870.3 

1469.5 

-1469.5 

4650.0 

1071.6 

841.9 

-841.9 

2100.0 

1850.6 

1454.0 

-1454.0 

4700.0 

1070.3 

841.0 

-841.0 

2150.0 

1829.8 

1437.7 

-1437.7 

4750.0 

1069.4 

840.2 

-840.2 

2200.0 

1808.2 

1420.7 

-1420.7 

4800.0 

1068.6 

839.6 

-839.6 

2250.0 

1785.5 

1402.9 

-1402.9 

4850.0 

1068.0 

839.1 

-839.1 

2300.0 

1761.7 

1384.2 

-1384.2 

4900.0 

1067.5 

838.8 

-838.8 

2350.0 
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Table  A2.  Dimensions  of  the  Corner-Fillet  Surfaces  of  the  Contraction  (Curve  2). 
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